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Abstract. Wc examine the radial asymptotic behavior of spherically symmetric 
Lcmaitrc-Tolman-Bondi dust models by looking at their covariant scalars along 
radial rays, which are spacelike geodesies parametrized by proper length i, 
orthogonal to the 4— velocity and to the orbits of SO (3). By introducing quasi- 
local scalars defined as integral functions along the rays, wc obtain a complete and 
covariant representation of the models, leading to an initial value parametrization 
in which all scalars can be given by scaling laws depending on two metric scale 
factors and two basic initial value functions. Considering regular "open" LTB 
models whose space slices allow for a diverging £, wc provide the conditions 
on the radial coordinate so that its asymptotic limit corresponds to the limit 
as £ — > oo. The "asymptotic state" is then defined as this limit, together 
with asymptotic series expansion around it, evaluated for all metric functions, 
covariant scalars (local and quasi-local) and their fluctuations. By looking at 
different sets of initial conditions, wc examine and classify the asymptotic states 
of parabolic, hyperbolic and open elliptic models admitting a symmetry center. 
We show that in the radial direction the models can be asymptotic to any one of 
the following spacctimcs: FLRW dust cosmologies with zero or negative spatial 
curvature, sections of Minkowski fiat space (including Milne's space), sections of 
the Schwarzschild-Kruskal manifold or self-similar dust solutions. 
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1. Introduction. 

The spherically symmetric LTB dust models pQ are among the best known and most 
useful exact solutions of Einstein'e equations. Since they allow us to examine non- 
linear effects analytically, or at least in a tractable way, there is an extensive literature 
(see [2J [3] for comprehensive reviews) using them, mostly as models of cosmological 
inhomogeneities [H |5j [H [3 [8], but also as in other theoretical contexts, such as 
gravitational collapse and censorship of singularities [lOj [TT] or quantum gravity [12] . 
There is also a widespread literature [I3J Ell ESI ESI EH EH] considering LTB models 
as tools to probe how the cosmic acceleration associated to recent observations can be 
accounted for inhomogeneities, without introducing an exotic source like dark energy. 
LTB models are also a standard choice [HI HOI HH to apply Buchert's scalar 
averaging formalism [23], in which the effects of dark energy could be mimicked by 
"back-reaction" terms (see [24 for a review of all this literature). In practically all 
articles the models are parametrized by a standard set of free functions and analytic 
solutions (implicit and parametric). 

While the literature is vast and exhaustive, there is still room to explore further 
development on their theoretical properties (see for example [25] in this context). The 
present article deals with a theoretical problem that has not been, as far as we are 
aware, previously examined in the literature, namely: the asymptotic behavior of LTB 
dust models (through their covariant scalars) in the radial spacelike direction, which 
can be defined in a covariant manner in terms of spacelike geodesies (radial rays) 
whose tangent vectors are orthogonal to the 4-velocity and to the orbits of SO (3). In 
the remaining of this section we explain and summarize the contents of the article. 

Basic background on LTB models is provided in section 3: the metric, field 
equations, classification in kinematic classes: parabolic, hyperbolic and elliptic, as 
well as a covariant time slicing that defines the space slices 3 T[t] orthogonal to the 
4-velocity field and marked by constant values of t. In section 3 we introduce an 
alternative set of quasi-local scalar variables [26l EHJ [27l [29l [30] , leading to a "fluid 
flow" description of the dynamics of the models that is similar (and equivalent) to the 
"1+3" approach of Ellis, Bruni and coworkers [3TJE2], which f° r LTB models (as with 
all spherically symmetric spacetimes) reduce to scalar equations |33j . As we showed in 
PTl [S!J1 [317] , these variables can also be understood in the framework of a perturbation 
formalism on a FLRW "background" defined by the quasi-local scalars, which satisfy 
FLRW dynamics, while their fluctuations are gauge invariant and covariant non-linear 
perturbations. The quasi-local variables and their fluctuations are potentially useful 
for a numeric approach to LTB models (see [55] and sections XI-XIII of [3D]), but they 
are also very handy for analytic and qualitative work, since they lead to an initial value 
parametrization of the models in which all covariant scalars can be given by simple 
scaling laws that depend on the two metric functions (scaled to an initial 3 T[^]) and 
initial value functions. 

In section 4 we discuss the definition and properties of the proper radial length, £, 
along the radial rays, which is the affine parameter of these geodesies. Since we need 
to explore how scalars behave as I — > oo, we will only consider LTB models (admitting 
a symmetry center) in which this limit can be realized, which means "open" models 
whose 3-dimensional space slices 3 T[t] orthogonal to the 4-velocity are homeomorphic 
(topologically equivalent) to K 3 , as £ is always finite at all 3 T[t] in "closed" elliptic 
models with slices 3 T[<] homeomorphic to S 3 . 

In order to probe the asymptotic behavior of covariant (local and quasi-local) 
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scalars along the radial rays, we should (ideally) evaluate these scalars as functions of 
t, which without numeric work is practically an impossible task (even qualitatively), 
as I is an integral of a metric function that can only be determined (at best) in implicit 
form in terms of another metric function. Instead, we provide the conditions so that 
this radial asymptotic regime can be examined in terms of the dependence of scalars 
on a well defined radial coordinate. Once this is done, the limit r — > oo will correspond 
to £ — ¥ oo, and as a consequence, a covariant characterization of the asymptotic regime 
can be provided in terms of the radial coordinate and initial value functions. 

Since we are considering LTB models that comply with basic regularity conditions 
(no shell crossings and scalars may only diverge at a central singularity), we define in 
section 5 a regular asymptotic regime in the radial direction by demanding that all 
scalars are smooth and finite as r — > oo (which now corresponds to £ — > oo). We also 
discuss (Lemmas 1 and 2) the relation between I and R — y/gee, which is another 
important invariant in spherically symmetric spacetimes [34] . Since the quasi-local 
scalars satisfy less complicated scaling laws than local scalars, it is more practical to 
study their radial asymptotics first. Hence, we prove in Lemma 3 in section 6 that 
both types of scalars, the local and quasi-local, share the same asymptotic behavior 
(given a common set of assumptions on the initial value functions). We also choose 
a radial coordinate gauge (as there is a radial coordinate gauge freedom in the LTB 
metric). 

In section 7 we consider the radial asymptotic behavior of the r-dependent initial 
value functions, which are the quasi-local scalars and their fluctuations evaluated at a 
fiducial "initial" slice 3 Ti — 3 T[<i] (the A qi and 8\ A ^), and are necessarily restricted by 
regularity conditions (the Hellaby-Lake conditions [HI OM 135] ) • We assume for these 
functions a uniform asymptotic convergence to specific (but not restrictive) asymptotic 
trial analytic functions of r: power law, logarithmic or exponential. However, before 
using these convergence forms to probe the scalars A q , A and the 8^ by means of the 
scaling laws derived in section 4, we introduce in section 8 the distinction between the 
"asymptotic limit", which is simply the limit of scalars as r — > oo, and the notion of 
"asymptotic state" , which we define as the set of these limits, together with suitable 
series expansions around them, evaluated for the metric functions and all covariant 
scalars and fluctuations. 

We examine the asymptotic limits and states separately for parabolic (section 9), 
hyperbolic (section 10), elliptic (section 11) models and in section 12 for special models 
with a simultaneous big bang (initial central singularity) and maximal expansion. A 
summary of these asymptotic limits and states, as well as a brief discussion, are 
provided in section 13. Depending on the initial value functions, the scalars in 
open LTB models have as asymptotic limit, either a FLRW cosmology (with zero 
or negative spatial curvature) or a section of Minkowski spacetime given in "non- 
standard" coordinates that generalize those defining Milne spacetime (the particular 
solution "[s2]" of [5]). Within those LTB models whose asymptotic limit is a Minkowski 
section we can recognize asymptotic states that clearly identify various particular cases 
of LTB models: sections of Minkowski spacetime that include the Milne universe, 
self-similar LTB solutions (see pages 344-345 of [3] and [TU1 [3S]) or sections of the 
Schwarzschild-Kruskal spacetime given in terms of coordinates constructed with radial 
geodesies (Lemaitre and Novikov coordinates, see page 332 of [3]). 

As we comment in section 13, LTB models asymptotic to a FLRW cosmology can 
be understood in the context of embedding a dust inhomogeneity in a dust FLRW 
background, without resorting to an artificial matching at a fixed comoving boundary. 
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In other words, these models can be fully relativistic and less artificial representations 
of the Newtonian "spherical collapse models" used as toy models of structure formation 
(see [29] for discussion on this point). This embedding into homogeneous cosmology is 
not compatible with LTB models asymptotic to Minkowski (in any of their asymptotic 
states). Instead, some of these models (specially the vacuum dominated hyperbolic 
ones) can be considered as toy models of a dust inhomogeneity surrounded by a 
large cosmic void. In this context, these configurations can approximate a spherical 
realization of the notion of "finite infinity" (or "fi") suggested by Ellis [37], and 
considered further by Wiltshire [55]) as a description of an intermediate scale in 
which galactic cluster structures can be studied as (approximately) asymptotically 
flat configurations. 

The article contains three appendices: Appendix A provides the analytic 
(parametric and implicit) solutions in the conventional variables, Appendix B discusses 
the possibility of considering local scalars (instead of quasi-local ones) as initial value 
functions, while Appendix C summarizes the various particular case spacetimes that 
follow from specializing the free parameters of LTB models (in our description: the 
initial value functions). As shown in sections 9-12, the radial asymptotic state of any 
open LTB model corresponds to a specific spacetime in this list. 



2. LTB models, kinematic classes and a fluid flow time slicing. 

LTB dust models in their conventional variables are characterized by the following 
metric and field equations G ab — Kpu a u b 

ds 2 = -c 2 dt 2 + dr 2 + R 2 (d6 2 + sin 2 6d<p 2 ), (1) 

1 + E 

R 2 = 2 ~f+E, (2) 

2M' = npR 2 R', (3) 

where p is the rest-mass density, n — 8irG/c 2 , E — E(r), M — M(r), while 
R = u a V a R = dR/d{ct) and R' = OR/ dr. 

It is common usage in the literature (see [5K3JIH]) to classify the solutions of ^ in 
"kinematic equivalence classes" given by the sign of E, which determines the existence 
of a zero of R 2 . Since E = E(r), the sign of this function can be, cither the same in 
the full range of r, in which case we have LTB models of a given kinematic class, or it 
can change sign in specific ranges of r, defining LTB models with regions of different 
kinematic class (see [SJ). These kinematic classes are 

E = 0, Parabolic models or regions (4a) 

E > 0, Hyperbolic models or regions (46) 

E < 0, Elliptic models or regions (4c) 

where the equal sign in (4b) and (|4c[) holds only in a symmetry center. The solutions 



of the Friedman-like field equation (|2| for each kinematic class are given in Appendix 
A. The case M = with E > arbitrary has been classified in [8] as the solution 
"[s2]" and is locally equivalent to Minkowski spacetime. 
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2.1. A fluid flow time slicing and covariant scalars. 

The normal geodesic 4-velocity in ([I]) defines a natural time slicing in which the space 
slices are the 3-dimensional Riemannian hypersurfaces 3 T[t], orthogonal to u a , with 
metric h ab — u a Ub + g a t, and marked by arbitrary constant values of t. Each spacelike 
slice is a warped product 3 T[t] — xM( r ) X R ^ 2 (0>V?)j where the warping function 
is R(t,r) > 0, the fibers are concentric 2-spheres S 2 with surface area AirR 2 (t : r) 
(orbits of SO (3)), while the leaves xM( r ) are "radial rays" or curves of the form 
C(r) = [cto,r,0Q, tpo], with to,9o,ipo constants. The rays are orthogonal to the fibers 
and isometric to each other in any given 3 T[i], and are also geodesies of the 3 T[t] 
and spacelike geodesies of the LTB metric Since we are assuming the existence 
of (at least) one symmetry center, then every xM( r ) f° r ^ constant is diffeomorphic 
to K + = {x | x > 0}. Hence, we will consider every LTB scalar function as equivalent, 
under the time slicing given by u a , to a one parameter family of real valued functions 
A[t] : R+ -> E so that A[t](r) = A(t,r). 

Besides p given by (|3|, other covariant objects associated with LTB models are 
the expansion scalar, 0, the Ricci scalar of the space slices, 3 1Z, plus the shear and 
electric Weyl tensors, a ab , E ab 

e -v M *-^ + ^ *n- 2{ERy (5) 

° ~ VaU - R n ~ R?R' ' (5) 

^ab = V (aMfo) - (Q/3)h ab = E ~ a \ E ab = u c u d C abcd = £ ~ ab , (6) 

where h ab = u a u b - g ab , V a = h b a V b , and C abcd is the Weyl tensor, E ab = h ab - 3n a r] b 
with rf = \J h rr 8r being the unit tangent vector along the radial rays (orthogonal to 
u a and to the orbits of SO(3)). The scalars £ and S in ^ are 



3 



R R' 
R~ Ri 



K 



Since LTB models (as all spherically symmetric spacetimes) are LRS (locally 
rotationally symmetric), they can be completely characterized by covariant scalars. 
Considering ([6| and ((7|), these are the local "fluid flow" scalars 

{p, 6, E, 8}, (8) 

whose evolution equations completely determine the dynamics of LTB models in the 
fluid flow or "1+3" approach [3TJ[32], and thus provide an alternative approach to 
that based on the analytic solutions of 



3. Quasi local scalars and their fluctuations. 

For every scalar function A in LTB models we define its quasi-local dual A q as the 
family of real valued functions A q [t] : R + —> R given by [|][f] 

= f AFdV p _ r AR*R'&x 
q f TdV p f'RZR'dx ' [ ' 

X This section provides the minimal background material to make this article as self-contained 
as possible. The reader is advised to consult reference [30| for details on the quasi-local scalar 
representation of LTB models. 

§ Since it is clear that t is a fixed arbitrary parameter we will omit henceforth the notation [t] unless 
it is needed. 
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where the integration is along arbirtary slices 3 T[i],with T = (1 + E) 1 / 2 , dV p = 
y/h a ijdrd8dip, and we are using the notation J Q r ...da; = f^~^ ...dx. The quasi-local 
scalars comply with the following properties 

4 = (^)' = ir i A -M> 



Mr) -Mr) 



1 



A' R 3 dx. 



qV ' R 3 (r)J Q 

Given the pair of scalars {A, A q }, we define the relative fluctuations as 

a 'j a i _ 1 r mta. 



§(A) = A M 

A„ 



3R'/R A q {r)R 3 (7 



A' R 3 dx, 



where we used (10a I and (106) 



(10a) 
(106) 

(11) 



3.1. Scaling laws and initial value functions. 

The quasi-local scalars lead in a natural manner to an initial value parametrization of 
LTB models, so that all quantities can be scaled in terms of their value at a fiducial (or 
"initial") slice 3 % = 3 T[ti], where t = t j is arbitrary. Hence, the subindex j will denote 
henceforth "initial value functions" , which will be understood to be scalar functions 
evaluated at t = ti. This procedure suggests rephrasing the metric functions R and 
R' as dimensionless scale factors 

£ = -, r=4^ = i + 4^, (12) 

leading to scaling laws for all covariant scalars, which now become functions of the scale 
factors L, T and initial value functions. Considering the definition ([9]), the quasi-local 
duals of p, and 3 1Z depend only on L and initial value functions 
m ql M 

M = — = (13) 



in,, 



~rJ 

kqi _ 

L 2 ~ 
a L 2 

1/2 _ 

q ~ L 2 



k, 



' i? 3 ' 
E 

^R 2 ' 

R 2 _ , _ 2m ql - 
R 2 ~ Zmq kq ~ L* 



kqiL 



(14) 
(15) 



where, to simplify the notation, we have introduced (and will use henceforth) the 
definitions 

« K , 3 K , 3 K q „, 6 

2m = — p. 2m„ = —£»„, k = , k„ = -, ft = — , 

3 q 3 Hq ' 6 ' q 6 ' 3 

Scaling laws for the local scalars follow readily from Q and ^ as: 



rt 



e„ 



(16) 



l 3 1 J l 3 r' 



fc=^ [! + *<*>] = ''' 



L 2 



L 2 r 



r - i 



3(1 + *< W ) 



(17a) 
(176) 



Since local fluid flow scalars in (|8|) are expressible in terms of m q , H q and k q and their 



fluctuations (11 1 
m = m q 1 + 5^" 



T~L — 1~L Q 



-m q 5 {r ' 



1 + 5M 



k — kn 



1 + 



(18a) 
(186) 
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we get, with the help from ( 11 ), ( 13 1, ( 14 1 and ( 15 ), the scaling laws for the fluctuations 

(19a) 



1 + 6^ = 



1 



- 4 m) 

r 



.*W = ?/3 + « 



2S m = 2m « 5(m) 



(196) 



kg <y(*> 



2m qi 5 {m ^ 



k qi LS^ 



2m 



2m qi [8, 



(m) 



2m„ 



i - r] - k qi l [6, 



(fe) 



Kg, Z/ 



(i-r)] 



which allow us to obtain scaling laws for the local expansion scalar, %, and the scalar 
function associated with the shear tensor, S. 

We note that the scalars m qi T-L q , k q and their fluctuations are covariant objects, 
as M, E, R, R = w a V a i? are invariants in spherically spacetimes [34]. Since the 
conventional variables M and E depend only on r, it is convenient to define them as 
initial value functions 



M = m ql R°, E 
Given Q and plj ), the LTB metric takes the form 

> 2 R? dr 2 



ds z 



-c z dt z + L 



1 kqi 



R 



sin 2 9d<tf 



(21) 



(22) 



The Friedman-like equation |2]) now takes the form (15). Its solutions, are equivalent 
to those of Q in Appendix A, but now expressed in terms of L, m q i, k q i. These 
solutions will be given explicitly in sections 9, 10 and 11. We can compute from these 
solutions a scaling law for T. 



3.2. Curvature singularities. 



The scaling laws ( 13 )— ( 20 1 clearly indicate the existence of two possible curvature 



singularities whose coordinate locus is 



L(t,r) — 0, central singularity 

T(t, r) = shell crossing singularity. 



(23a) 
(236) 



so that for reasonable initial value functions (bounded and continuous), all scalars 
A q = m q , k q , T-L q diverge as L — > 0, whereas local scalars A = m, k, H can also 
diverge if T — > (even if L > 0). Notice that if T > 0, then all scalars A and A q 
only diverge at the central singularity L = 0, which is an intrinsice feature of LTB 
models. However, if T — > for L > 0, then all the relative fluctuations 5^ diverge 
(with A q ^ 0) , so that local scalars A diverge while their quasi-local duals A q remain 
bounded. This is an obviously unphysical effect of shell crossings that must be avoided. 
We will denote by "regular LTB models" all configurations for which shell crossing 
singularities are absent, thus complying with 

F > V (ct, r) such that L > (24) 

In order to test this regularity condition we need to compute T, which will be done 
separately for parabolic, hyperbolic and elliptic models in sections 9, 10 and 11. 
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4. A well behaved proper radial length. 

In order to examine the radial dependence of scalar functions we need to define the 
"radial direction" in precise and covariant terms. There is no inherent covariant 
meaning in the radial coordinate. In fact, the metrics ([I]), ( |22[ ) and h a b are invariant 
under arbitrary re-scalings r = r(f), indicating the existence of a coordinate gauge 
freedom that can always be used, either to simplify computations or to eliminate any 
initial value function by using it as radial coordinate. Yet, radial rays are geodesic 
curves whose affine parameter is radial proper length, and so the radial dependence 
of scalars at every individual 3 T[t] acquires a covariant meaning by relating it to this 
parameter. 

The proper radial length along an arbitrary 3 T[t] can be defined as the function 
£[t] : R+ -> R such that 

t[t](r)= / y/g^dx = / — dx = LT-^dx, (25) 
Jo Jo •> Jo •> 

where T = Vl + E = [1 - k ql Rj] 1/2 , (26) 

so that ^[i](0) = for all t. A well behaved proper length must necessarily satisfy 

(r) > and l'{r) > for r > at all 3 T[i], so that r 2 > n <s> £(r 2 ) > £{r x ) at all 



3 T[t]. Considering (23a) and (24), these requirements are satisfied for regular radial 
rays in regular LTB models if the following regularity condition among initial value 
function holds 

sign {Rl l (r) ) = signer ). (27) 
Thus, if Ri has a zero, it must be a common zero of T (and of the same order). While 
proper length is the natural parameter to characterize radial dependence in a covariant 
manner, it is not convenient to use it as a spacetime coordinate because in general: 
£ = £(t,r) (the same remark applies to another important scalar like R). Hence, for 
practical reasons we need to describe the radial dependence of scalars in terms of their 
radial coordinate. Since we can always use radial coordinate gauge freedom to fix Ri, 



(27) can be understood as a consistency condition on the radial coordinate, so that it 
effectively mirrors the dependence on proper radial length along radial rays. 

Since (25) is valid for all 3 T[i], then for regular LTB models (L > 0) the sign of 



R' is the same as the sign of R[ for all r. Therefore, as long as (27) holds, the radial 
coordinate can be used to probe the radial profiles of scalars at an arbitrary 3 T[i], as 
these profiles will be qualitatively analogous to those with respect to £: 
f)T? f) f? /?' 

R ' = ~dl £ ' = ^l^ ^ S ign(R')= S ign(dR/d£), (28a) 
r)A f)A W 

A = ~M^' = ~My sign(A')-sign(cM/d£). (286) 

While a zero of A' necessarily corresponds to a zero of dA/d£ at any individual 3 T[t], 
in time dependent scalars A = A(t, r), a zero of A' can arise (in general) at different 
values of r (or £) for different 3 T[t], or it can arise in some of the 3 T[i] and not in 



others, all of this without violating (24), (27) or (28a|— (286). As we mentioned above, 



the case of R is different: a zero of R' must be common to a zero of J 7 , and so it 
necessarily occurs at a fixed r and is common to all 3 T . 



Notice that (27) could be violated even if shell crossings are absent ( (|24j) holds). 
In these situations there would be a surface layer at the comoving worldline marking 
the zero of R\ [HI US] ■ It is evident that such a surface layer also implies an ill-defined 
proper radial length. 
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5. A regular asymptotic regime in the radial direction. 

A radial asymptotic regime is associated with convergence and behavior of metric 
functions and covariant scalars the limit I — > oo along radial rays in the space slices 
3 T(t). However, given the existence of (at least) one symmetry center, these slices 
can be homeomorphic to either R 3 ("open" models) or to S 3 ("closed" models), which 



must be elliptic because of (27 1. Since I is everywhere finite in closed models, then 
the limit £ — > oo can only be realized in open models, which can be of all kinematic 
classes (hyperbolic, parabolic or elliptic with — 1 < E < Q). We will only consider 
open models for the remaining of this article. 



Assuming regular open models for which (24 1 and (27) hold, we will consider 
the radial asymptotic regime as regular if the scalars {m, H, k} and {m q , 1-L q , k q } 
converge to finite values as £ — ¥ oo. If we adopt as a criterion to define a curvature 
singularity that curvature scalars diverge if a given coordinate locus is reached by 
geodesies in finite affine parameter values, and since £ is such a parameter, then there 
would be no singularity (technically speaking) if these scalars diverge as I — > oo. Also, 
it is possible to conceive a situation in which k, H and their quasi-local duals diverge 
(even at finite £) while the densities m, m q remain bounded, so that a singularity does 
not arise because curvature scalars remain bounded (see [8] and Appendix A of [30]). 
While these situations cannot be ruled out, we will exclude them from consideration 
and will, nevertheless, assume henceforth that all covariant scalars (to, fc, % and their 



quasi-local duals) may only diverge at a central singularity (23a I, remaining bounded 
everywhere and also in the limit £ — > oo. [JJ 

5.1. Relation between £ and R. 

Since R' > holds everywhere in open models and R is an important invariant scalar 
involved in the definition of £, it is necessary to examine the relation between the 
limits £ — > oo and R — ¥ oo . We prove now the following 



Lemma 1 . Let £(r) given by (25) be the proper radial length along an 
arbitrary 3 T[t] homeomorphic to M 3 in a regular LTB model, the limit 

lim R{£) = oo, (29) 

l— foo 

holds for all parabolic and hyperbolic models, and for elliptic models in which 
J- = \J\ + E converges to a nonzero constant as £ — > oo. 



Proof . As a consequence of (|28a|), we have 

~d£ 



which is valid at each 3 T[i] separately. Since T = 1 for parabolic models and 
T > 1 and T' > hold for all regular hyperbolic models, the result follows 



directly from (30). For open elliptic models, T lies in the range < T < 1 



with J"(0) = 1. If P < then as £ -» oo we have T -> J" with < F Q < 1, 



then (30) implies 

Tq I < R{£) < £, (31) 

|| The relative fluctuations <5'' 4 ' might diverge for finite r under regular conditions, see Appendix A4 

of ED] 
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and so the result follows. If J- — > 1 as I — > oo, then T' must have a zero at 
some r = r tv (which is a minimum of J 7 because J-(0) = 1), then T > ^"(rtv) 
holds for all r ^ r tv and J 7 ' > for all r > r tv - In this case constraining 
bounds similar to (31) can be constructed from ( |30[ ) with J-q = J-(r tv ) and 
the result follows. However, if T — >• as I — > oo, then (depending on how T 
converges to zero) R might converge to a finite constant in this limit. 



Corollary. If the limit (29 1 holds in one 3 T[t], it must hold in all 3 T[t] 
proof is straightforward, since Lemma 1 is valid for arbitrary 3 T[t]. 



The 
Also, 



the 3 T[t] constitute a smooth foliation of LTB models and (if standard reg- 
ularity holds) the rays in all 3 T[i] are complete geodesies in the direction of 
increasing £. 



Lemma 2: (converse of Lemma 1). The limit 

lim £(R) = oo (32) 

R— ¥00 

holds for all parabolic and elliptic models, and for all hyperbolic models in 



which k„ 



constant as R — > oo. 



Proof . Since T is related to k q by (14) and (21), we can rewrite (25) (and 
also (30)) along an arbitrary 3 T[i] as 



<(R) =/V 



dR 



k q R 2 ]V* 



(33) 



where we now consider k q = k q (R) and we used the fact that dR = R'dx is an 



exact differential because t is kept constant as the integral in ( 25 1 is evaluated. 

In parabolic models we have k q = everywhere, so the result follows trivially 

from (33). For open elliptic models, (33) implies that k q R 2 < 1 must hold 

1/2 > i 



for all R, which is equivalent to l/J 7 
have in general 

t(R) > I dR = R 



k q R z 



Therefore, we 



and so 



if R — > oo. For hyperbolic models, (33) takes the form 



[l + |fc ? |i? 2 ]V2' 



(34) 



(35) 



It is evident then, that the integral in (35) converges only if \k q (R)\ diverges 



hence £(R) will always diverge as R — > oo if \k q (R)\ — > constant in this limit. 
Notice that this result is valid at individual (but arbitrary) 3 T[£]. In general, 
\k q (R)\ would converge to a different constant in different 3 T[i]. 



6. Radial asymptotics of local and quasi local scalars. 

Bearing in mind that we are only considering open models for which R' > holds 
everywhere, then a convenient (and practical) way to choose the radial coordinate is 

R i (r) = R Q r, (36) 



where Rq is a constant length scale. As a consequence of (36), radial dependence 



becomes dependence on the initial value function Ri and Rq provides a characteristic 
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length scale for the radial coordinate. Since (36) complies with (27), then r — > oo 
implies I — > oo in all open models (including elliptic models for which R tends to a 
finite constnt in this limit) . Also, as a consequence of the corollary of Lemma 1 , r — > oo 
also implies R — ¥ oo for all open models (save the above mentioned elliptic models). 
Unless specified otherwise, we will henceforth assume that the radial coordinate has 



been fixed by the gauge ( 36 1 



Assuming the coordinate gauge (36) and that (24) and (27) hold, we examine now 
the relation between the radial asymptotic behavior of local and quasi-local scalars. 

Lemma 3 . The following result holds in any given 3 T[t] 

lim A{r) = A lim AJr) = Aq, (37) 

r— s-oo r— >-oo 

where Ao is a constant. 

Proof . We consider the case when A' > for sufficiently large values of r. 
The case when A' < is analogous. If the limit of A as r —> oo is Aq, then 
for all e > there exists u(e) such that Aq — e < A(r) < Aq holds for all 
r > u(e). Constraining A by means of this inequality in the definition of A q 
in ([9]) leads immediately to Aq — e < A q (r) < Aq, hence A q — > A as r — > 00. 



The converse result follows from (10 a 



The following results follow trivially from Lemma 3: 
Corollary 1: A' q — > and A' — > both hold as r — > c 
Corollary 2: A q S^ = A - A q ^ a,s r ^ 00. 



Since we have assumed that (27) holds, these results are valid as I — > 00 along 
individual but arbitrary 3 T[i], though the constant A will be different at different 
3 T[i\ and (in general) A = A (t). These results are useful because it is easier to probe 
the asymptotics on the quasi-local scalars A q first, as they satisfy less complicated 
scaling laws (that do not involve T) and the analytic solutions are given in terms of 
m q i and k q i. Once we have worked out the A q , the asymptotic behavior of the A 
follows from Lemma 3 and its corollaries. 



6.1. The relative fluctuations 8^ A \ 

Lemma 3 yields the following general result: 

Corollary 3 of Lemma 3: If A — >• Ao 7^ (or, equivalently, A q — > Aq ^ 0) 
as r — > 00 in an arbitrary 3 T[i], then 5^ — > in this limit, irrespective of 
how fast or slow A and A q converge to Aq. The proof follows directly from 
Corollary 2 of Lemma 3. 

However, if A — ¥ as r — ¥ 00, then this corollary no longer applies. The limit value 
of 5^ in this case is not (necessarily) zero, but depends on the specific asymptotic 
form in which A and A q decay to zero. 



7. Asymptotics of initial value functions. 

The results proven so far are valid at individual (but arbitrary) 3 T[t]. However, it is 
very difficult to actually test them because we need to evaluate the involved scalars 
as functions of r for arbitrary fixed t. With the exception of parabolic models, this is 
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very hard because L is, in general, given implicitly in the solutions of the Friedman- 



like equation ( |15| (equivalent to ([2])). However, the coordinate gauge (36) provides a 
simple relation between Ri and r, and this allows us to obtain general analytic results 
on the asymptotic behavior of the initial value functions Ai, A q i and their fluctuations 
8\ A \ As we show in sections 9, 10, 11 and 12, these analytic expressions yield analytic 
forms for the asymptotic behavior of these variables for t^ti. 

7.1. Asymptotics of the 5\ A ^ 

In order to obtain the limit of the when A, A q — > Aq = 0, we need to make specific 
assumptions on the convergence of these scalars to zero. 

Uniform asymptotic convergence. Let Ai and Ai be smooth integrable scalar 
functions on 3 T[i;] both tending to zero as r — > 00. The scalar A4 uniformly 
converges asymptotically to Ai if for every e > there exists a real positive 
number y, such that \Ai — A\\ < e for all r > y. We denote this convergence 
by Ai ~ Ai. Notice that this definition is also applicable to the quasi-local 
scalars A q i converging asymptotically to a function A q i. 

We can assume a given asymptotic convergence either for Ai or for A q i. We 
show bellow how prescribing A qi we obtain Ai and &\ A \ The alternative approach 

(prescribing Ai to obtain A q i and 5\ A ^) is discussed in Appendix B. 

It follows directly from the definition of uniform asymptotic co nver gence tha t if 
A qi ~ A qi then A' qi ~ A' qi holds for r > y. Therefore, considering (10a) and ( Jllj ), if 
we assume A qi ~ A qi we obtain 



A * ~ Aqi + = A ** + ^r 1 > ( 38ffl ) 



A'qi rAL 
ZR'JRi 
A A) A' qi /A qi _ f A' qi 



3R'JRi 3 A 



(386) 



where we eliminated Ri from ( 36 ) . 



7.2. Initial density 

While a non-negative m q i is an initial value function of LTB models of all kinematic 
classes (parabolic, hyperbolic or elliptic), the sign and regularity conditions associated 
with the quasi-local spatial curvature, k g i, depends on the kinematic class. Hence, we 
examine in this subsection the admissible forms of asymptotic convergence for m q i, 
leaving the discussion for the convergence of k q i for sections 10 and 11 that deal with 
the asymptotics of hyperbolic and open elliptic models. 

For whatever form we might choose for the asymptotic convergence of mi or m q i, 
we must be very careful that both of these densities remain non-negative, which (if 



standard regularity holds and from the scaling laws (13) and (17al) implies that m q 
and m are also non-negative at all 3 T[i\. Since m q i appears as initial value function 
in the analytic solutions in our parametrization of the models (see [30]), we examine 
specific asymptotic convergence forms for it. 
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Considering (38a I, the asymptotic convergence forms for m q i 
and m. 



in which both 



remain non-negative everywhere are: 
Logarithmic: m q i ~ mo/ \n r a , (a > 0) M 
Power law: m g j ~ m r~ a (0 < a < 3), M 

0, we have m q 

m 



For the power law (396 1 wit h a 

mo and 6^ 
3 also implies m, - 



0. For q39&p with a > and for ([39a 



m i?i]r 3 /lnr Q , (39a) 
m Rlr 3 - a . (396) 

mo? and so Lemma 3 yields 
>, q i — > 0, and so Lemma 



0, but the asymptotic behavior of rrii is not identical. We have 



from ( 38a H 386 1 



Power law: 



Logarithmic: 



to 



('-!) 



a 
3, 

too (3 In r 
3a In 2 



1) 



(m) 



f(m) 



a 

3' 

1 



31nr 



(40 a) 

(406) 
becomes 



Notice that a power law with a > 3 in (396) implies that rrii in ( |40c 
negative for sufficiently large r, though it remains close to zero and (from Lemma 
3) it still tends to zero as r — > oo. The same would happen for an exponential decay 
m qi ~ TO exp(— ar). This undesirable behavior is consistent with the asymptotic 



forms for M in (39a)-(396), since standard regularity requires this function to be 



monotonously increasing, thus admitting for m q i only a power law decay with a < 3 
or a slow logarithmic decay. The fact that rrii < happens for steeper decays of m q i 
follows because rrii ^ m qi holds for r > y if the radial gradients are negative in this 
range, thus, if m q i becomes very close to zero in a very steep decay it may force rrii 
to become negative (but close to zero). 

While not all assumptions on the decay of m q i to zero yield a positive local density 
rrii ; ah assumptions on the decay of rrii — > necessarily yield rrii and m q i positive in 
the full range r > y. This is so because m q i > rrii > holds in this range (as and 



m' qi are negative, see (10al). We discuss this issue in Appendix B. 



8. Asymptotic limits and asymptotic states. 

We explore now the asymptotic behavior of regular LTB models along radial rays of 
arbitrary 3 T[i] (i.e. t =^ t{ arbitrary and finite), and looking at parabolic, hyperbolic 
and elliptic models separately. This involves evaluating asymptotic series expansions 
of the metric functions and covariant scalars around their limit r — > oo, under 
the assumption that the initial value functions m q i and k q i converge to admissible 
forms of uniform asymptotic convergence summarized in tables 1 and 2, and given 



explicitly by (39a)-(396) and (see section 10 and 11) (59a)-(596| and (94a)-(946) 



Since every LTB model can be completely determined by m q i and k q i (assuming 



the gauge (36)), once we assume admissible convergence forms for these initial value 



functions, we can find asymptotic convergence forms and asymptotic series for the 
metric functions L, T, the scalars m q , T-Lq, k q , their fluctuations S^ m \ 5^ , 5( k \ as 
well as local scalars to, H, k and other auxiliary quantities (M, E, t bb , etc). All this 
information completely characterizes an "asymptotical state" for every class of models 
based on a specific assumption on the convergence forms of m qi and k q i and the 
kinematic class (parabolic, hyperbolic and elliptic). 

The resulting limits and expansions in the characteristic functions contained in the 
asymptotic states will be compared to the equivalent parameters in those spacetimes, 
listed in Appendix C, that are particular and limiting cases of LTB models, namely: 
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dust FLRW cosmologies, Schwarzschild-Kruskal, Minkowski (including Milne) and 
self-similar LTB solutions. The radial asymptotics of specific classes of LTB models 
can be then characterized on the basis of this comparison, which we express below 
more precise terms: 



Let {M.,g} be an LTB model manifold with metric g given by (22), with 
{M. (o) , 5(o) } being the spacetime manifold and metric of any one of the par- 
ticular cases listed in Appendix C. Let X = {L, T, A, A q }, where A, A q 
denote the scalars m, k, % and m q , k q , % q in {A4,g}, with Xf \ being the 
equivalent set of functions in {A^(o),.9(o)}- 

Definition. We shall say that {M, g} is radially asymptotic to {A^(o), 3(o)}i or 
that {A4,g} converges asymptotically to {M(o)i9(o)} m the radial direction, 
if the series expansion for every function in X around its limit as r — > oo 
coincides with its equivalent function in A( j up to leading terms. 

Since it is clear that we are considering asymptotic behavior along radial rays, we 
will simply state that a given LTB model or class of models is "asymptotic to" or 
"converges to" a given particular case. Notice that it is sufficient to evaluate the limit 
as r — ¥ oo for the quasi-local scalars A q , as (from Lemma 3) this limit will be the 
same for the local scalars A. 

However, the evaluation of the strict limit as r — > oo of the functions in X is not 
sufficient to characterize the asymptotic behavior of a given class of LTB models, since 
the same limit can correspond to different asymptotic expansions and convergences. 



For example, the limit as r — > oo of the metric ( 22 ) of some LTB models can coincide 
with a Minkowski metric in spherical coordinates (L,T — > 1), but if we consider 
the expansions of L and T around this limit, the metric of these models could become 
equivalent (up to leading terms in the expansions) to the metric of an LTB self-similar 
solution or a section of Schwarzschild-Kruskal spacetime (whose metrics themselves 
have a Minkowski limit as r — > oo). In order to distinguish these cases, we need to 
consider all metric functions and scalars in X, evaluate their asymptotic expansions 
(asymptotic state) and compare them (up to leading terms) with the functions and 
X( ). Notice that Another point that needs explaining is why we did not include 
the fluctuations in the set X in the convergence criterion above. The reason is that 
A — > (or equivalently A q — > 0) as r — > implies (in general) a nonzero limit for 
5^ A \ while the equivalent of 5^ in A( ) could be identically zero. This can happen 
when A(q) is associated with a Minkowskian particular case, for which m — m q = 
holds identically, but then we have m q and m —¥ in an LTB model converging 
to this particular case, and thus does not tend to zero. However, the asymptotic 
convergence is well characterized because the nonzero limit of 8^ is consistent with 
m q — ¥ as r — ¥ oo. A summary of the asymptotic behavior for different classes of 
LTB models obtained in sections 9-12 is given in section 13. 



9. Radial asymptotics of parabolic models. 



For k q i — the Friedmann-like equation ( 15 ) yields a closed analytic expression for L 



L = 



1 + -^/2m ql c(t - U) 



2/3 



(41) 
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which is equivalent to the parametric solution (A.l), and we are only considering 



expanding configurations (L increases for t > tj). The bang time, t bh , follows by 



setting L = and t = t bb in (41 1 

ct hh = ct 



3y/2m qi 



Cti 



By differentiating (41) and from (12 1, we obtain 



1 



i-(rri) 



L 3 / 2 ' 



(42) 



(43) 



Since k a 



are m q i and S. 



ki — and Ri has been fixed by ( 36 1 , the only initial value functions 
^ Necessary and sufficient conditions to fulfill ( 24 1 (absence of shell 



crossings) are given by the Hellaby-Lake conditions, which for parabolic models and 
in terms of our initial value functions takes the form |30j 

(44) 



1 < si m) < 0, 



which from ( 11 ) implies m' qi < 0, so that m q i tends to zero or to a nonzero constant and 
its admissible asymptotic forms (preventing a negative TOj) are given by (39a)-(396). 



Assuming an asymptotic convergence m q i ~ rh q i given by either one of (39 a 



(396) the asymptotic convergence for L follows from (41 1 as 

o 12/3 



L~ L = 



1+2 V 2 ™qi c ( i _ *») 



Considering (13), (15), (21 1 and (42 1 we get 



^27 



' qi 



IF' 



(45) 
(46 a) 

These asymptotic forms are valid for any regular parabolic model, while (from Lemma 
3) the limits of m, H as r — >■ oo are the same as those of m q , T-L q . Also, Lemmas 1 
and 2 imply that R — > oo for all regular initial value functions. The asymptotic limits 
and states depend on the choice of m qi . Figures la and lb depict the domain in the 
(ct, r) plane parabolic models with asymptotic limit to FLRW and Minkowski. 



M ~ M = fh q i Rq r 3 



Cti - 



9.1. Parabolic models asymptotic to spatially flat FLRW. 



If 



r(m) 



oqi —r too > 0, so that m q i = Too (the power law form (396) with a = 0), then 



as r — > oo follows from the corollary 3 of Lemma 1. We have the asymptotic 



convergence forms 



L(t) = 



-i 2/3 



u i~ u « ~ L s/2 (t y 

M - M = too Rl r 3 , 
so that 



l + -y/2rr*c(t-ti) 



r - i 



to 



ct v , 



L 3 {t)' 

Ctuu = Cti 



6 (m ) - 0, 



3V2too 



= const. 



L -> L(t), f4l, m q ->m q (t), H q ^H q (t), 6™ -> 0, -»• 0. 



(47a) 
(476) 
(47c) 

(48) 
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Since 8^ and 8^ 



tend to zero and F — > 1, then to m„ and H. — > "H„ (which 



is consistent with Lemma 3). It is evident that the asymptotic state from these 
expansions is consistent with the corresponding parameters of a spatially fiat dust 



FLRW model (see equations (C.1)-(C.3|) 



9.2. Minkowski limit and its asymptotic states 
If TO 



as r — > oo, then (45) and (43|) yield 
L~L«1 + y/2m q i c(t - U), r~f«l 



3- 



~8 { i ' n) ^M^ l c{t-t i ). 



(49) 



where we have written 
(1 + e) 2 / 3 » 1 + (2/3)e, 



(instead of 
valid for e <C 1. 



: ~") because we used the approximation 
The power law decay of m q i is useful to 
illustrate the asymptotic convergence forms for the remaining quantities. Considering 

ri q ~ 



with < a < 3, we obtain 

^2TO^r-"/ 2 



ft 



8. 



(m) 



Z 3 / 2 

-a/3, 



M ~ M = to ffi r 



toq r 



L 3 



: -a/3, 
3V2too 



— > — oo, 



a 
3' 



(50 a) 
(506) 

(50 c) 



where we used (19a) and (20). From Lemma 3 we have to — > and H — > 0, which 
agrees with (18a). If we only consider the limit as r — > oo, then TO g i — > clearly implies 
a Minkowski limit, since L — > 1 and T — > 1 transforms the LTB metric ( p2| (with 
fc g i = 0) into a Minkowski metric in spherical coor dinates. The fact that 8^\ 8^ 
tend to a nonzero value is consistent with to 9 , 1-l q tending to zero (as corollary 3 of 
Lemma 3 does not apply). If m q i decays logarithmically as in (39a), we obtain the 



same Minkowski limit, but with fluctuations 8^ A ' now also tending to zero. 
The asymptotic states of parabolic LTB models with fh q i — > are: 

• Asymptotic to spatially flat self similar solution. If a = 2, then L = £(C), where 



( = c(t — ti)/r is the self-similar variable in (C.17). It is evident that the metric 
and all functions asymptotically converge to their respective forms in the spatially 



flat self similar solution given by ( |C.14[ ), ( |C.15[ ), ( |C.17| )-( |C.19| wit h k = 0, which 
is the case E = in equation (2.29) of [36 (see also pages 344-345 of [3] and 
[TO]). Notice that parabolic LTB models with this asymptotic behavior are not 
self-similar solutions, they only converge to a self-similar solution as r — > oo 
along the 3 T[t]: the latter solution only follows if m q i — too r~ 2 holds exactly in 
all the domain of r. 

Asymptotic to Schwarzschild-Kruskal. 



If 



3, then a comparison with 



(C.10), (C.ll) and (C.12) implies that we have an asymptotic convergence to 



Schwarzschild-Kruskal solution in spatially flat comoving coordinates (Lemaitre 
coordinates, see page 332 of [3]). This convergence is consistent with M ~ 
TO i?o = constant, as opposed to M diverging for a < 3. Notice that 8^ — s- —1, 
hence to = m g (l + 8^) converges to zero much faster than m q ~ m i?o/i? 3 . 

Asymptotic to Minkowski in generalized Milne coordinates. If a / 2, 3, then 



comparison of L in (49) and (C.6) shows an asymptotic state compatible with 



that of a locally Minkowski spacetime in coordinates that generalize Milne's. 
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However, the form (49) does not not follow from the particular solution with 

M = m„ 



' qi 



in (A.6). Rather, we have an asymptotic state compatible with a 
section of Minkowski, but one based on expansions around an asymptotic limit 
m q i = in models (parabolic) for which m q i is not strictly zero but k q i = 
holds everywhere. This Minkowsky asymptotic state occurs also for hyperbolic 
and elliptic models converging to a parabolic model (see sections 10 and 11). 



10. Radial asymptotics of hyperbolic models. 



Analytic solutions for hyperbolic models follow from solving the Friedman-like 
equation ( 15 ) for k q i < 0. These solutions are equivalent to those of ([2| given by 
(A.2) and (A.2) in parametric and implicit form. Wc will work with the implicit form 
given in terms of our variables by 

(f> = Zh(xiL), 

where is the function 

u i ^ Z h {u) = u 1/2 (2 + u) 1/2 - arccosh(l + u). 

and 

<j) =(j}{t,x i ,y i )=y i c{t-t i ) + Z h {x i ), 

\kqi\ l&qil^ 2 



(51) 

(52) 

(53a) 
(536) 



The bang time emerges from setting L — in (51) as the following function of m q i 
and I Awl: 



ct bb — cti 



Zhjxj) 
Vi 



(54) 



The metric function T follows from (|51|) by implicit derivation since it is related to 

3 



L'/L by dT2l). The result is 



r = l + 3(^ m) -6l K >) ( 1- 



rlq 

Ha 



-m q c(t-ti) [8, 



:(m) 



(55) 



where "H 9 and H q i follow from (151, while c(t — U) is given by (51) and (53a). The 
Hellaby-Lake conditions [5J O 35 to fulfill the condition ( 24 ) for absence of shell 
crossings are given in terms of initial value functions by |30| 

4 ^ (m) > -i, 



where 



3RyRi 



H 



c(U - t bb ) 4 



qi 



3 A( fc ) 



(56) 
(57) 



follows by differentiating (54). 
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rhqi 


IM 


Xi 


Section 10.3 
Xi -» 


Section 10.4 

Xj — > oo 


Section 10.5 

Xj -> x 


PL 


PL 


x Q r a -$ 


< a < P < 2 


< /3 < a < 3 


< a = /3 < 2 


PL 


LOG 


Xq r a 1 lnr' 9 


a = 
P > 


< a < 3 
/3 > 




LOG 


PL 


xq r~° lnr a 


< P < 2 
a > 


/3 = 
a > 




LOG 


LOG 


Xq a//3 






a,/3 > 



Table 1. Combinations of the admissible asymptotic convergence forms 
for hyperbolic models. The table presents the asymptotic convergence of 
%i = \ kaiWf hqj for the four possible combinations of the forms ih q i and \k q i\ in 
l |39a| |-( [39fe| and ( |59a| l-| |596| ). The terms "PL" and " LOG " stand for "pow er law " 
and "logarithmic" decay forms, respectively given by ( |39a| l, |59a| and ( |39fe[ l, |59fe|. 
The behavior of Sj examined in sections 10.3-10.5 follows from the restrictions 
on the parameters a, j3, hence the results of these sections can be readily applied 
to initial value functions having each asymptotic convergence form. 



10.1. Initial negative curvature 



It is evident from (11), (21) that the profile and asymptotic convergence of the initial 
value function k q i < in a hyperbolic model is restricted by the second regularity 
condition in ( 56 1 , which for a negative k qi necessarily implies that E = J- 2 — 1 = 
—kqiR^r 2 must be monotonously increasing (and that fcj < 0). In order to examine 



the case k q i < 0, we remark (see (11)) that 
sign(c^ fc >) = sig]i(*i/fc 9 ). 



(58) 

Hence, 5^ has the same sign in the combination k q < 0, k' q > as in k q > 0, k' q < 0. 



< holds everywhere (with possibly ^ ql 







As a consequence, if we assume that k, 
asymptotically), then we can examine negative spatial curvature by applying (40a)- 
(406) directly to \k q i\. Considering all these points together with (21) and (36), the 



asymptotic forms of \k q i\ and E compatible with standard regularity conditions and 



with a well defined asymptotic radial range (r — > oo =>■ £j 



1 </' 



fco/ln r, 

k a r~P (0<p<2), 



0. 



E 



fc Rl r 2 



oo) are 
/ In r, 



3' 



E - k R, 



2,H 



(59a) 
(596) 



where fco > is a constant. Notice that \k q i\ ~ ko (or j3 — 2) allows for both an 
increasing k' qi > or decreasing k' qi < behavior, though (from Lemma 3) k' qi — > as 



r — > oo and \k q i\ — > fco- In both cases (59a) and (596) we have \k q i\ decaying to zero. 
The asymptotic forms of the local curvature |fcj| are analogous to those of \k q i\ given 
above, since the allowed decay forms are not steeper than r -3 (see Appendix B). 
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10.2. Asymptotic approximations for the analytic solutions 



There is no closed explicit form for L as (41 ), hence we must work with the implicit 



solutions (51), with the correspondence rule of the function Zh given by ( 52 ) . Since the 
admissible forms in (59a)-(59o) imply that \k q i\ is bounded in the radial asymptotic 



regime, then the radial asymptotic behavior of <f> = Zh(xiL) depends on the behavior 
of Xi in this regime, which in turn depends on the choice of initial value functions m q j 
and \k q i\. Since Zj t is a monotonously increasing function of its arguments, then (51 1 
implies 







Xi — » OO Zh(Xi) — > oo =>■ (p — > 

while, another possibility is given by 

Xi — >• xq — const., 4> {kqil 1 ^ 2 xq c(t 







Z h (x,iL) - 
Z h {xi~L) 



0. 



U) + Z h (x Q ) Z h (x L). 



(60 a) 
(606) 

(61) 



The three patterns outlined by (60a)-(606) and (61 ), and listed in the 4th, 5th and 6th 



columns of table [T] provide all possible combinations of asymptotic behavior of regular 
hyperbolic models in terms of their initial conditions (through x{). It is possible to 
apply appropriate approximations for each one of these cases in order to obtain L, 
which follows by inverting (|51|) 



L 



(62) 



Given an asymptotic form for L obtained in terms of admissible convergence forms of 



m q i and \k q i\ in (39a)-(39o) and (59a)-(59&), we can obtain asymptotic convergence 



(or approximated) forms for the remaining quantities. We examine separately the 



asymptotic behavior associated with each of the three cases (60a)-(60&l and (61), 
each one listed in the fourth, fifth and sixth columns of table [TJ Figures la and 
lb respectively depict the domain in the (ct, r) plane of hyperbolic models with 
asymptotic limit to FLRW and Minkowski. 



1 0.3. Matter dominated asymptotics: convergence to parabolic models. 



If m q i dominates over | k q i \ 



in the limit r — > oo, then x-, — > and we have the situation 



described by (60a). Hence, we can consider the following approximation for Zh 

/2 



Zh(u) 



u 3/2 _ V* u 5/2 



for 



0. 



(63) 



3 20 

where we use the sign to distinguish this approximation from the approximations 
that follow from assumptions on uniform asymptotic convergence ( "~" ) on the initial 
value functions. Applying (63) (up to the leading term) to Zh(xi) and Zh(xiL) in 



(51), and considering that m 



i ql and \k qi \ 



]^2rh qi c(t - U) 



\k q i\, yields 

2/3 



(64) 



which coincides with the form of L in ( 45 ) , though now L does not (necessarily) 



converge uniformly to L, but approximates it up to a certain order related to the 
convergence of Xi to zero (hence this approximation depends also on assumptions on 



k q i\, which is strictly zero in parabolic models). Nevertheless, (64) indicates that 



hyperbolic models with initial conditions in which Xi — > converge to a parabolic 
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model. Though, we have now \k q i\ ^ 0, hence the asymptotic state might not define 
a general parabolic model and the forms of the involved quantities should not be 
identical to those of parabolic models in the previous section. In order to discuss 
these issues, we examine below the two types of asymptotic states associated with 



( 64 1 that emerge from the form of rh q 



Asymptotic to spatially flat FLRW. 

If m q i ~ mo > (the power law form (396) with a — 0, see top entry in the 
fourth column of table B, then S^ m) -> and Jiil yields L » L{t) with L{t) 



given by (47a), which is the scale factor of a spatially flat FLRW dust model. 
The conventional variables M, E and t bh have the asymptotic forms 

2 



M ~ m R^r 



3„.3 



E 



\k qi \R r 2 , 



Cti 



3 v / 2m^ 



const. 



(65) 



(59a 



where \k q i \ can take the power law form (59 6 1 with /3 > or the logarithmic decay 
Since L{t) is constant along the 3 T[<] and \k qi \ — > 0, then the asymptotic 



form s of th e re m ainin g time -d epe ndent scalars follow from series expansions of 



(131, (14), (151, (19a), (1961, (20) and (551 around \k qi \ ~ \k qi \, leading up to 



order 



k q i\ tO 



III,, 



r 



m 
L 3 : 



\k. 



L 2 



"►0, H a 



1 



"> 1, 



28™ 



3(L 5 / 2 -l)|fc gi |ff fc 
4m i 3 / 2 

(2 + 3% fc) )(L 5 / 2 -l 
4m L 3 / 2 
3(L 5 / 2 -l)\ L\k. 



y/2m~Q 

3(L^~l)\k ql \o^ 
4m £ 3 / 2 



y/2m 
Z 3 / 2 ' 



(66 a) 



ki\~4 k) 



0, (666) 



(66c) 




(AO 



2i 5 /2 



0. 



m 



(66 d) 



where \k q i\, 6^ can take any of the admissible asymptotic forms (|59a|)-(|596|) 



The fact that, as r — > oo, we have m q , T-L q tending to time dependent forms, 
while r — > 1 and 5^ tends to a perturbative term, clearly shows that these 
hyperbolic models are asymptotic to the spatially flat FLRW dust cosmology. 

Asymptotic limit to Minkowski. 

k q i\ — > 0, following either power law decays with < a < j3 or 



If rh q i — > and 

any combination of power law and logarithmic decay, as given by (39 a 
p^-([59l, we 



the case PL-PL with a 
( 64 ) becomes as r — >• oo 

L 



J39&p and 

have the cases listed in the fourth column of table [1] (excluding 



0) . For all these combinations of m qi —> and | k 







1 + y/2m qi c{t-U) -> 1, 

where we used (1 + e) 2 / 3 « 1 + 2/3e, which holds for e <C 1 
variables M, E and t bb take the asymptotic forms 

2 



M 



|A: 9l |i?o?" 2 , 



(67) 

The conventional 
(68) 



Zy/2rh qi 

Hyperbolic models with these initial value functions converge to parabolic models 
whose asymptotic limit is Minkowski (see previous section), but they do not 
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converge to a general parabolic model of this type. This point can be b etter 
illustrated by considering a power law decay for both m q i and \k q i\ as in (396) 
and ( |596[ ) (see the case PL-PL in the top entry of fourth column of table [l] with 
a > 0). Since we must have a < f3 for Xi — > to hold, then there cannot be a 
convergence to parabolic models for which 2 < a < 3 holds because regularity 
conditions imply /3 < 2. 



The scalars m q and \k q \ and T-L q take the forms 



o. 



L 2 



o, M 

m„ 



tin 



^ ( 1+ tJ v 91 



1 - ^2rh q ic(t - U) + ~ 



0, 



(69 a) 



(696) 



where L is given by (67). Considering the expansions above and only terms linear 
v„i or Xi, and since yl + Hi/2 w 1 + holds for Xj <C 1, we obtain up 



m ^/TOq 

to leading order 



r w 1 + 5y2rf^c(t - tj) ->■ 1, 



?(m) 



f« » *<« - iy^-c(< - t,) tfj ro) (2 + 3^ (fc) ) -+ *f \ 



(70 a) 
(706) 
(70c) 



25W « # m) - ^^2^c(t - ^} m) (l + s| ro) ) - - ^ (fe) )^ -> ^ (m) . (70d) 



where the asymptotic convergence form of <5>j m \ depends on the choice of 



m gl , \k qi \ in (39a|-(396) and (59o|-(596l 



10.4- Vacuum dominated asymptotics. 



Consider now the case when \k q i\ dominates over m qi in the limit r 



oo, so that 

Xi — > oo and we have the situation described by (606), with initial value functions 



listed in the fifth column of table [T] We can consider then the following approximation 
for Z h 



Zh(u) 



ln(w) 



for ti>l, 



(71) 



where, as in the previous subsection, we distinguish this approximation from that 
following a given assumption of uniform asymptotic convergence on the initial value 
functions (hence the sign "«"). Applying (71) (up to the leading term) to Zh(xi) and 
Zh(xiL) in (51), and considering that m q i 



ih q i and \k ql \ 



\k qi \, yields 



Lw 1 + M 1/2 c(t-ti), 



(72) 



so that L approximates the solution of the equation L 2 = which is the Friedman 
equation (15) with m q = and k q < 0, or equivalently, the solution of ([2| in (A.6) 
with M = and E > 0. From (IC?5j) (fcl9j) in Appendix C2, these LTB models 



are the solutions [s2] in [8], and are locally Minkowskian (equivalent to sections of 
Minkowski spacetime parametrized by non-standard coordinates that generalize the 
Milne universe) . Evidently, hyperbolic models characterized by initial value functions 
complying with Xi — > oo converge in the radial direction to these Minkowskian sections. 
Notice that k q < does not imply that spacetime curvature is nonzero, as k q is the 
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quasi-local dual to the spatial curvature (fe) of the 3-dimensional slices 3 T[i\. The 
3-dimensional scalar curvature of these hypersurfaces can be non-trivial, even in a 
vacuum flat Minkowski space. On the other hand, m q is associated with a source m 
appearing in T ab and generating spacetime curvature. If \k q i\ dominates over m q i in 
the asymptotic radial range, then the frame dependent kinematic effects of the spatial 
curvature play the major role in the behavior of all incumbent scalars in this range. 

The spatial curvature dominated hyperbolic models lead to the following 
asymptotic states depending on the convergence of the initial value function \k q i\: 



Asymptotic to a Milne Universe 
while m 



qi\ ko > (which corresponds to the power law form (59 b I with 



ql —r has any of the forms ( 39a I— ( 396 ) with a > 0, then (72) becomes 
L(t) = 1 



L PS L(t) = 1 + y/k c{t-ti), 

while the conventional variables take the asymptotic form 



0). 

is 

(73) 



M — rh qi R r 



E ~ k R 2 r 2 , 



cti- 



1 

/ko 



const. 



(74) 



By comparing with (C.7)-(C.9), it is evident that hyperbolic models with these 
initial value functions converge to the locally Minkowskian Milne universe. 

Bearing in mind that L = L(t) and m q i — » 0, then the asymptotic forms of 
the remaining time-dependent scalars follow from series expansions of (13), (14), 
(196), (20) and (55) around rh qi = 0. We obtain up to leading order 



(151, (19a 



in m™ the following forms simi 



o. 



ko_ 



ar to (|66a|)-(|66d) 



1 



fen L 



v/fcp 

L : 



3c(t — U)d- in 



(™) ^,2 
qi 



2&W 



■3/2 



?(m) 



3c(t 



u)~sl m) (1 



Urn) 



■,3/2 



[i + Vhc{t - u)f 



2c(t 



x( m ) ™2 
qi 



U)S- Hi 



■.3/2 



[1 + Vh>c(t - u)} 
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(™)^ 



in a 



k [1 + Vkoc(t - ti)] 



(75 a) 



(756) 



(75 c) 



(75 d) 



(75 e) 



;(m) 



can take any of the admissible asymptotic forms (39 a)- (39 6) 



where rh q i, 
with a > 0. 

Asymptotic to a generalized Milne universe. 

If xi diverges asr-> oo, but both rhqi, \^qi\ given by (39a)-(396) and (59a)-(596) 
tend to zero in this limit, we have the case listed in the fifth column of table |l| (of 
course, excluding the Milne universe discussed above, hence a, f3 > necessarily 
holds). As a consequence, L takes the form (72), which implies L ~ L — > 1 as 
r — > oo. The functions M and E take the forms in (68) with the appropriate 
values of a, (3 > 0, while i bb takes the asymptotic form 

ct bb ps ci, - — -> -oo, (76) 
\kqi\ 1 



Radial asymptotics of LTB dust models. 



23 



Hyperbolic models with these initial value functions converge to generalized 
versions of the Milne universe (c ase [ s2] in [5]). The scalars m Q a nd k q are given 
by the same expressions as in (69a I, but with L given by (72) and rh q i, \k q i\ 
complying with ii = \k q i\/rh q i oo. The scalar H q in (15) takes the following 
asymptotic form 

1/2 



IM 1/2 



i 



Xi L 



IM 1/2 



(l + Xi) 
->0, 



Ifc^l 1 / 2 [l-\~k qi \ 1/2 c(t-t l )+x l 
1 1 

Considering only terms linear in Im^l 1 / 2 and Xi, we obtain up to leading order 



ith 



(77a) 
(776) 



l+~c(t-t i W>\% 9i \ 1 '* 



1. 



6 {k) 



5 



(m) 



1/2 



?(m) 



c(fe) 



c(t- tl )(2 + 3f)f Ifc,,! 1 /^^, 



2$"° -5< fc) )x, 



(78 a) 
(786) 
(78c) 
(78d) 



where <5>j , 5^ depend on the choice of m qi , \k qi \ in (39a)-(396) and (59a)-(596). 
See figure lb. 



• Asymptotic to Schwarzschild-Kruskal. If we choose the power law decay for rh q i 
in (396) with a — 3, then 

M <~ tuoRq = const., 



m 



91 



UIqT 



(79) 



then a comparison with (C.10), (C.ll) and (C.12) reveals an asymptotic 
convergence to Schwarzschild-Kruskal solution in coordinates given by geodesic 
observers with positive binding energy (see page 332 of [3]). Evidently, the 
hypersurfaces 3 T[i] of the Schwarzschild-Kruskal spacetime itself, in this specific 
time slicing, has an asymptotic limit to Minkowski in the ra dial d i rectio n. The 
asymptotic form for \k q i\ can be any of the admissible forms (59a|-(596), but in 
order to illustrate the asymptotic behavior of the scalars we will also assume a 
power law decay \k q i \ ~ \k q i\ = fco with (3 < 2, so that 

^--l + Ofe), ^--f + Ofe), (80a) 

r /3/2 

ct hb cti = — > — oo. (806) 



moa-3 



k 

E ~ k Rlr 2 ~P 



< \k a 



The asymptotic forms for the remaining scalars simply follow from specializing 
the parameters of (77a)-(78c) to this case: 



rig 

r ? 

S (r> 



ykp 

„/3/2 



1 - 



m 



rfS/2 



il _ ^c(<-,,) +0(r _ a 
« -i + 0(r -fl ) -> -1, 



fc r 3 -/ 3 



(81a) 

(816) 
(81c) 
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6 
3 



r P/2 

2to 



•0(r" 




3' 



k r 3 ~P 3 ' 



24 
(81 d) 
(81 e) 



Notice that (81c) implies that local density, to, decays to zero much faster than 
7Ti g , as expected in a convergence to Schwarzschild-like vacuum conditions. 

.Z0.5. Generic hyperbolic asymptotics. 



We consider now the case (61) when x% — > xo = constant, listed in the sixth column 
of table [l] We will assume the power law decay (396) and (596) with a = fi, so that 
density and spatial curvature decay at the same rate (the treatment of the logarithmic 
decay is analogous). This decay implies m q i ~ too r~ 7 and \k q i\ ~ fcor -7 , where 
< 7 < 2 (regularity rules out 7 > 2). Hence, we obtain from (51) 

fcp c(t - tj) 

where y = y knXn and is given by (52). Thus, following (51) and (52), we have 
L ~ L so that ( 62 ) becomes 

2/0 € + Z h (x ) ~ Z/,(afo-£(0)i ( 83a ) 
1 



cc 



(82) 



x 



Z^{y t + Z h {x )). 



(836) 



where we can write L = L(£) because Z^{xq) is a constant. We have then for 7 > 



lim L = lim L(0 = — Z^(Z h (x Q )) = 1. 



(84) 

As a consequence of ( |82| , (83a)-(836) and (84), and bearing in mind that we are 
considering t constant, the asymptotic form of L for very large r is follows from a 
series expansion of Z(£) for £ <C 1 up to order £ (i.e. order 0(r~ 7 / 2 )): 

1(0 « 1(0) + L i£ (0) £ = 1 + v^mo + fc 6 (85) 
where we used (84) and i.j(O) = [dL/d£]j =0 was computed by implicit derivation of 
fl83a| ). 

The following asymptotic convergence forms follow from (13), (14), (15), (21) and 

fro 



(54): 



H 



[2m + fc L} 1/2 



too 



r7 /2 L 3/2 y r7i 3' 

M - TO ^r 3 " 7 , E ~ k Rlr 2 -~<, ct b 



ct 



r iL 2 ' 

Z h (x )r^ 2 
Vo 



(86 a) 
(866) 



Since S^ m) , <*f° -> -7/3 and considering ([ssj, then Jl9a| ), |l96|), (|55| and (666} lead 
to the following approximations up to order r -7 '/ 2 

(87a) 
(876) 
(87c) 
(87 d) 



r re 1 


- ^ (2to 



+ M 1/2 e, 




7 
3 


1 - 


3 6 7 (2TO + fc ) 1/2 £ 


! 




7 
3 


1 - 


2 ; 7 (2TO + fc ) 1/2 d 

6 


5 




T 


1 - 


7[2TO (3- 7 ) + fco(2 


-7)] 


3 


6 (2777 + k ) V2 
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The remaining scalars m, k, % are approximately equal to m q , k qi H q times a factor 
1 



7/3. If instead of a power law form we consider a logarithmic as in (39 a I and 

£j)/lnr and <5 Z 



(59a I, we obtain the same forms, except that £ = c(t 



*F-> 0. 



The asymptotic expressions (85), (866|-(86a) and (87a)-(87d) hold for any regular 
hyperbolic model for which X4 = \k q i\/m q i — > Xo = constant. Since the free parameter 
is 7, we can identify the following asymptotic states: 

• Asymptotic to negatively curved FLRW. If 7 = 0, then we have m q i ~ mo and 
\k q i\ ~ fc , as well as £ = c(i — ij). Hence, L, % q and m g depend only on t, while 



asymptotically Mar and E oc r and £ bb — > constant. Comparing with (C.l), 



(C.2) and (C.4), we can see that L, M, E and t bb converge to their corresponding 



forms of a FLRW dust model with negative spatial curvature. 

Asymptotic to generalized Milne. If < 7 < 2 then £ —> as r 

—00 and r — > 1 



?7lq, |fcq 



and Tig tend to zero, while t h 



00, thus 
These limits and 



the expansions around it, together with (C.6) clearly point out that models with 



these parameters are asymptotic to a section of Minkowski in generalized Milne 
coordinates. 

Asymptotic to the self-similar solution with negative spatial curvature. If 7 = 2, 



then m q i, \k qi \ oc r~ z and all quantities in (85), (866)-(86a) and (87a)-(87d) 
converge to their respective forms in the self-similar solution with negative spatial 
curvature in (C.14) (this is the case E > in equation (2.29) of [36 ). 



11. Radial asymptotics of open elliptic models. 



As with the hyperbolic models, we will use the implicit analytic solution of ( 15 ) for 
k q i > 0, equivalent to (A.5), which in terms of our variables takes the form 

Z e (xiL) expanding phase 

(88) 

27r — Z e (xiL) collapsing phase 

where the expanding and collapsing phases correspond to H q > and H q < 0, and Z e 
is given by 



with 



u i-> Z e (u) = arccos(l - it) - u 1/2 (2 - u) 1/2 



cj> = <f>(t,Xi,yi) = yi c(t - U) + Z e (xi), 



A- 



y, 



3/2 

qi 



(89) 

(90 a) 
(906) 
The times 



iii q i 1 'veil 

The scale factor L is restricted by < L < L max = = 2m gl /fc gi . 

associated with the initial singularity, the maximal expansion (L max ) and the collapsing 
singularity are given by 

Z e ( x i) K 27T 

c£ bb — cti , ci max = ci bb -\- , c£ co n — ct hh -\- . (91) 

Vi Vi Vi 

Notice that (in general) i max = t max (r) and t coll = t coU (r), so (like t bh (r)) neither one is 

(in general) simultaneous (see section 12 for the cases with simultaneous t bb or i max ). 

For every comoving observer r = const., the time evolution is contained in the range 

t bb (r) < t < t coll (r). 
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The metric function F takes the same form as (55), with T-L q given by (15) with 



01301133] to comply with (24) are 



k qi > and c(t — t{) follows from (88) and (90a). The Hellaby-Lake conditions 



"&R[/ Rj 



<0, 



3i?^ / R, 



> 0, 



where ct' roU follows from differentiating ci coll in (91) 



ct' 



SR'i/Ri 



Am) 



?0) 



c(^coii ^bb) H~ 



> -i, 



bb 

3it^/ i?i ' 



while ct' bb has the same form as (57) with ct bb given by (91 ) 



(92) 



(93) 



11.1. Initial positive curvature. 

The asymptotic forms for k q i compatible with standard regularity must comply with 
the constraint 1 + E = T 2 = 1 — k q iR\r 2 , and thus k q ir 2 cannot diverge as r — > oo. 



Mi 



Hence, the admissible asymptotic forms are 



k ql = k a r- p 08 > 2), 



I 1 ~ 1 - fcni?, 



2^2-/3 



3' 



-/3r 



T 2 - 1 - k Q Rlr 2 e 



-0r 



S 



(AO 



(94 a) 
(946) 



form for r — > oo, hence fc a 



ho 



Notice that the case k q i ~ kg (i.e. (94a) with ft = 0) is not admissible as an asymptotic 



ds in this limit for all open elliptic models (though 



m 



— > mo = const, is allowed). We also remark that the logarithmic decay is also 



ruled out because k qi must decay at least as fast as r . Also, from (176) and (196), 



if 5 t - < — 1 then we have k < even if k q > in all 3 T[t]. This situation occurs if 



> 3 in (94a) and in (946) for which S, 



:(k) 



-oo (see Appendix B of [30 ) . 



11.2. Asymptotic approximations for the analytic solutions. 

As with hyperbolic models, there is no closed explicit form for L and we must work 



with the implicit solutions (88)-(91), with Z e given by (89). As opposed to the 



hyperbolic case, we have now two branches, but <f> = Z e (xiL) in (88) is a monotonously 



increasing/decresing function of XiL in the expanding/colapsing branches, so it has a 
branched inverse 



iL{<t>) = 



(95) 



However, from (89), (22) and (36) for open elliptic models with k q i > we must have 
<Xi<2 



and 



1 - k ql r 2 R 2 > 



k qi r 2 < l/R\ 



(96) 

so that Xi is bounded and k qi r 2 cannot be an increasing function as r — > oo (which 



is consistent with (94a)-(946)). Hence (as opposed to the hyperbolic case) <f> must 
remain bounded. 

As a direct consequence of | 94o[ )-(946l and (96), the asymptotic convergence 



fco = constant is not allowed. Therefore, open elliptic LTB models are 
incompatible with a positively curved FLRW asymptotic limit. This is not surprising, 
since these FLRW dust models are necessarily closed (the 3 T[t] homeomorphic to § 3 ). 

As with the hyperbolic models, the asymptotic behavior of open elliptic models 
depends on the initial value functions through Xi. If we take into consideration the 
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m ql 


h 


Xi 


Section 11.3 
Xi -» 


Section 11.4 
ii — ?> xo < 2 


Section 11.5 

Xi — > Xq = 2 


PL 


PL 




(3> a 
< a < 3 


2 < /3 = a < 3 


2 < /3 = a < 3 
fco = 2toq 


PL 


EXP 


Xq r Q e _/3r 


p > 
< a < 3 






LOG 


PL 


Xo r~P In r a 


/3 > 2 
a > 






LOG 


EXP 


xq e _/3r lnr a 


/3 > 
a > 







Table 2. Combinations of the admissible asymptotic convergence forms 
for open elliptic models. The table presents the asymptotic convergence of 
Xi = k q i/r h q j for the four possible combinations of the forms rh q i and k q i in 
( |39a,[ l ( |396| t and | |94a| -| |94fe[ l. The terms "PL" and " EXP " stand for "pow er law " 
and "exponential" decay forms, respectively given by ( |396[ l, l |94a| and | |39a[ l, |946|. 
The behavior of Xi examined in sections 11.3-11.5 follows from the restrictions 
on the parameters a, /3, hence the results of these sections can be readily applied 
to initial value functions having each asymptotic convergence form. 



restrictions (94a)-(946) and (96), then (88) and (95) allow for the following limiting 
regimes 



0. 



xq < 2, 



(97) 
The 



which will be discussed in detail, respectively, in sections 11.3 to 11.5. 
combinations of the admissible asymptotic convergence forms compatible with ( 96 ) 



that correspond to each of these cases are listed in 4th, 5th and 6th columns of table 
[2] For reasons that will be explained in the following subsections, the collapsing branch 
will be needed only in the case Xi — > 2 in (97) (see figure 2). 



11.3. Convergence to parabolic models. 



qt to dominate m qi in the limit r — > 00, though 



Since Xi < 2 there is no possibility for k, 
the opposite situation is possible and leads to the limit Xi — > corresponding to the 
first one of the cases in ( 97 1 and listed in the fourth column of table [2] As we show 
further ahead, the collapsing phase of (88) will not be needed (see figure 2), while Z e 
is monotonously increasing in the expanding phase. We have then 







Z e {xi) -> 







Z e ( Xi L) ->■ 0. 



Hence, we can approximate Z e in (89) by 

V2 



Z e (u) 



u 3/2 + V2 m5/2 



for 



3 20 

which, up to the leading term, coincides with (63) 
hyperbolic case, we assume m q i 



0. 



(98) 



(99) 



i q i and k q i 



qt ■ 



As in the den sity dominated 
and apply (99) to Z e {xi) and 
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Z e (xiL) in the expanding branch of (88), leading exactly to the same form for L 
given by (64 1 

2/3 



L 



1 + -^"V c (t - *f) 



(100) 



which indicates that open elliptic models with initial conditions in which x% — > are 
asymptotic to a parabolic model. 

While open elliptic and hyperbolic models with Xi —¥ have the same radial 
asymptotic behavior, their time evolution is radically different, as it is restricted by 



for all choices in (94a)-(946), the maximal 



the initial and collapsing singularity (see figure 2). Conside ring ( 99 ) and the fact that 
we are assuming xt — ¥ and k q 
expansion and collapse times ci ma 



ct,. 



= Cti 



Ctj 



x , ct co n given by (91) take the asymptotic forms 



k q { X i 

2tt - Z e (xj 
k q { Xi 



Cti 



Cti 



V% -1/2 
v2 _l/a 

■-T Xi 



k 



V2 - 



Qi 



2tt 



k qi Xi 



— > 00. 



(101a) 



(1016) 



where ih q i and k qi take the admissible forms compatible with Xj, — > 0. Since the locus 
of maximal expansion marked by ci max (r) corresponds to a maximum L = £ max = 2/xj, 
we have L max — >■ oo as r — ¥ oo and so: i? max = RorL max — > oo in this limit. 

Since the collapsing phase corresponds to ct > ci max (r), the fact that ct — ct max (r) 
diverges as r — > oo implies that the radial asymptotic range for all hypersurfaces 3 T[t) 
occurs within the expanding phase, which justifies the fact that we did not need to 
consider the collapsing phase of ( 88 1 for the study of the radial asymptotics in this 



case. This is clearly illustrated by figures 2a and 2b. Since open elliptic models have 
the same asymptotic behavior as in hyperbolic models in which Xi — > 0, the same 
asymptotic limits arise: 

• Asymptotic to spatially flat FLRW. 



If m q i ~ toq > (the power law form (39&I with a = 0, see top entry in the fourth 
column of table 2), then (100) yields L w L(t) with L(t) given by (47a), which is 



the scale factor of a spatially flat FLRW dust model. The conventional variables 
M, E and t bb have the asymptotic forms as in (65) (with — k qi instead of 
while the maximal expansion and collapse times, ci max and ct coll , are given by 



(101a) and (1016) under the specialization (396) with a = 0, so that Xi — > 0. The 



asymptotic forms of the remaining time-dependent scalars are readily computed 



as in the hyperbolic case. The result is exactly the same forms as in (66a)-(66d). 
Asymptotic to Minkowski. 

If ifiqi — ¥ and k q i — > 0, following either power law decays with < a < ft or any 
combination of power law and exponenetial decay, as given by (39a)-(|396[) and 



( 94a )-( 946), we have the cases listed in the fourth column of table 2 (excluding 
the case PL-PL with a = which corresponds to the spatially flat FLRW case). 
For all admissible combinations of initial value functions ( 100 ) becomes as r — > oo 



L w 1 + \j2m q i c(t 
while the conventional variables M, 



-U) -> 1, 
E and ct h 



(102) 

take the same asymptotic forms 



as (68), while c£ max and ci coll follow from ( 101a)-( 1016), and the scalars m„ and 
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Xi, we obtain up to leading order 

u, * ^— ( i - - 1 « y2^~ 



fcq have the same forms as in (69a). Considering only terms linear in 



T 



1 - ^2rh qi c(t - U) 



Xi 

4 



->0, 



(103) 



The remaining quantities have the sa me a symptotic forms as (70a)-(70d), but 



with kqi and d 4 now given by (94a)-(946). 



11.4- Generic elliptic asymptotics. 

We examine now the second case in ( |104| | listed in the fifth column of table [2j 

1/2 

Xi — > x n = const. < 2, 4> — > k qi xo c(t — ti) + Z e (xo) — > Z e (x L), (104) 
From table 2, the only combination co mpat ible w ith X j —> x$ is that in which both, 



m„ 



m„ 



and kqi, have the power law forms (39 6 1 and (94a) with same exponent 



mor 



k r 



2 < 7 < 3, x 



< 2. 



m 



(105) 



where the restriction on 7 follows from the possible common exponents in (39b) and 
(94a). Considering that Z e {xo) < it for all xo < 2 (we examine the case xq = 2 in 



the following subsection), the maximal expansion and collapse times for these forms 
of rh q i and k q i take the following asymptotic forms 



rf,. 



Cti 



cU 



7T - Z e (x ) ^ 7/2 
27T - ^ e (xo) 



OO, 



7/2 



(106a) 



(1066) 



As a consequence, only the expanding phase in ( 88 ) is needed to examine the radial 
asymptotics in this case (just as with the case Xi — > 0). If Xi — > xq, then as r — > 00 the 
scale factor L at maximal expansion tends to a finite value: L = L max = 2/xt — > 2/xq, 
though we have i? max = RorL^^ — > (2Rq/xq) r — > 00 in this limit. Since 2/xq > 1 for 
Xo < 2 and L = 1 at the initial hypersurface 3 T[i;], then this hypersurface necessarily 
lies in the expanding phase: U < i max (see figure 2b). 



Since we are only considering the e xpa nding pha s e in (88), the equation 



Z e (xoL) has exactly the same form as in (82) and (83a)-(836), with Z e instead of Zh 
Thus, p4| holds and we have L 
takes now the form 



1 + \J 2m - k £, 



1 as r — ¥ 00, while the asymptotic expansion (85) 
c{t-U) 



r7 /2 



(107) 



The asymptotic convergence forms for M, E, t h 



k q and Ti. q follow from (13), 



(15), (21) and (91), and are very similar to those in (|866|-(|86a 

-k Rlr 2 -\ 



M 



moRlr^, E 
[2m - k L\ 1 / 2 

r7 /2 ^3/2 



ct h 



Cti - 



Z e (x ) 



0, m q 



mo 



0, kq 



where yo = vfco^o and L is given by (107). Since 6, 



> —7/3 and considering 

(107), then ( lgal^J 196), (55) and (666) lead to the same approximations (up to order 



yo 

kp 
riE 2 
) S W 



0. 



(108 a) 



(1086) 



' 2 ) as (87a)-(87d|, but with k Q replaced by — fc . Since 7 = is not possible, then 



all open elliptic models complying with ( 105 1 have a Minkowski asymptotic limit. The 
following asymptotic states emerge: 
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Asymptotic to the self-similar solution with positive spatial curvature. If 7 = 2, 
then rriqi, k q i oc r~ 2 and M oc r and E ~ — fco-R 2 , hold. Hence, there is asymptotic 
convergence to the self-similar solution with positive spatial curvature in (C.14) 
with self-similar variable £ = c(t — ti)/r = ( (this is the case E < in equation 
(2.29) of [35]). 

Asymptotic to Schwarzschild— Kruskal. If a = 3, then M ~ tuoRq — constant. 



Comparison with (C.lOl, (C.ll ) and (C.12) reveals an asymptotic convergence to 



Schwarzschild-Kruskal solution in coordinates given by geodesic observers with 
negative binding energy (Novikov coordinates, see page 332 of [3J). 



11.5. The case xq = 2. 

The particular case Xi ~ xq = 2 is characterized by the same initial value functions in 
(105) with k — 2m (sixth column of table [2]). Equation (88) can be rewritten as 

c(t - U) 



Z e (2L(0)~7r±2v^e, £ = 



7-7/2 



(109) 



where the + and — signs respectively correspond to the expanding and collapsing 
phase, and we used the fact that Z e (2) = ir. 

There is a fundamental difference (in comparison with the case Xq < 2) in the 
locus of the maximal expansion time, ci max . Considering that both m qi and k qi decay 
as r then we can write in general Xi ~ 2 + 0(r~ 7 ). Hence, we have for x% 2 



Z e { Xl ) rj 2v / 2r _7/2 



1 



l+^r^ /2 + 0(r- 7 ) 



which, from (91), leads to 



Cth 



Cti 



Cti 



| [7T-Z e (x t )]r 7 / 2 
2\ 2m„ 
[2^-Z e (x,)] 
2v^7^ 

2\/2^ 



1 



V TO o 



(110) 



(111a) 



— t 



-00, 



(1116) 
(111c) 



so that t hb and t 
asymptotic line of t 



tend to curves that are symmetric with respect to the constant 
(see figure 2c). 



Therefore, as a consequence of (111a 



(111c 



there exist an extended collapsing region for open elliptic models for which 
x = 2. As opposed to the cases X{ — > and Xj — > x < 2, there is now a radial 
asymptotic range for the hypersurfaces 3 T[t] in both the expanding phase and the 
collapsing phase ci max < ct < ct coa . In fact, models in which xo = 2 are asymptotic to 
the models with a simultaneous i max that will be examined in the following section. 

If xi — > 2, then: L = L max = 2/xj — > 1 as r — > 00, but i? max = R^rL^^ R^r = 
Ri — > 00 in this limit. Notice that (111a) implies that the locus of the maximal 
expansion itself lies in the expanding phase. The asymptotic limit of L follows by 



inverting ( 109 1 



2Z~Z e - 1 (7r±2v / 2~W£). 



(112) 

1 in this limit (for 



Since £ — > as r — > 00 and n = Z e (2), then ( 112 ) implies L ~ L 
both the expanding and collapsing phases) . While this is the same limit as in the case 
xo < 2, the asymptotic expansions like ( |107 ) and the equivalents of (87a|-(87q) must 
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be computed up to £ 2 , since now \/2mo — ko — 0. The result (valid for the expanding 
and collapsing phases) is 



i(o « i(o) + \i ,me = i - 



(113) 



where L^(0) = [d 2 L/d^ 2 ]^ = o was computed by implicit derivation of (109) with 



respect to £. The asymptotic form for T-L q follows readily from (108 6 1 with ko = 2m, 



and (113) 



ro £ 



TOqC 



1 



3TOo 



0. (114) 

= at t = ti, 



r7 /2 13/2 l-3m £ 2 /4 
The asymptotic form for T cannot be computed from ( 55 ) because £ 
hence H q i « 0. Instead, we use the definition of V in (12) applied to L in (113) 

7 m £ 2 /2 



1 



Bearing in mind that 5, 
readily 

26^ 



r_U_ 

L 

(m) _ 



1 



1 - 



i 



(7 



l)m ^ 2 



1. 



(115) 



m £ 2 /2 2 
= —7/3, the following asymptotic forms follow 



7 
3 



1 

3 



(3 - 7) toq 2 
2 5 
(2 - 7) to 2 
2 € " 
(2 - 7) to 



1 



-0(£ 3 ) 
0(£ 3 ) 
£ 2 + 0(£ 3 



_7 
3' 

7 
3' 



27 
3 ' 



(116a) 
(1166) 
(116c) 



Hence, the case xo — 2 yields the same asymptotic states as xq < 2, but with a faster 
decay 0(£ 2 ) to the Minkowski limit. The models are also asymptotic to Minkowski in 
generalized Milne coordinates when 2 < 7 < 3, with 7 = 2,3 corresponding to models 
asymptotic to self-similar and Schwarzschild-Kruskal spacetimes. 



12. Simultaneous big bang or maximal expansion. 

So far we have assumed that the initial curvature singularity (L — 0) is not 
simultaneous, but given by the curve [ci bb (r),r] in the (ct,r) plane. While a constant 
t bb is incomplatible with parabolic models, non-trivial and perfectly regular hyperbolic 
and elliptic models exist for which t' bb = (see [71 [H [26]). As shown in 30], the 
Hellaby-Lake conditions to avoid shell crossings are simply (56) and (92) with t' bb = 0. 
Also, the locus [c£ max (r),r] of the maximal expansion time (% q — 0) in elliptic models 
or regions is, in general, not simultaneous, though regular models exist with ct' max = 
[7J 130) . We examine in this section the radial asymptotic behavior of regular models 
with these characteristics. 



12.1. Simultaneous big bang. 

Let t bb = t(o) denote the constant time value associated with L — 0, then from (54) and 



(91), the initial value functions m qi and k qi are necessarily linked by the constraints 

Zh(\k ql \/m qi 



c(U - i( )) = F h {m qi , \k qi \) = m qi 
c{U - t( )) = F h (m qi , k qi ) = m q , 



IM 3/2 



Z e (k qi /m qi ) 

7,3/2 
K qi 



hyperbolic, 
elliptic, 



(117a) 
(1176) 
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where the functions Zh and Z e are given by (52) and (89), and we have used (53a) and 
(90a) to express Xi and yi in terms of m q i and k qi . We can prescribe one of the pair 
and find the other by solving the constraints (117a) and (1176). While these 



rr> 



qi 7 k q i 



constraints cannot be solved analytically, we examine below the radial asymptotics of 
models with t bb = t(o) by looking at them qualitatively. 



Since the functions m q i, k q i that solve ( 117a) and (1176) for given values c(i,-— i(o)) 



define the level curves of Fh(m q i,\k q i\), we need to verify first if the admissible 
assumptions of asymptotic convergence of m q i , k q i are compatible with these functions 
lying in a level curve, and then if their limit as r — > oo (in the (m q i, k q i) plane) fixes the 
value of c(U — i(o))- Considering these assumptions, the following mutually exclusive 
cases emerge: 



m q i ~ mo, k q i r*> ko. These assumptions are only possible in hyperbolic models, 
hence we have Fh(m q i, \k qi \) ~ Fh,(ko,mo) — constant, which fixes the level curve 
c(ti — t(o))- Hence, all hyperbolic models with these initial value functions and 
^bb = t(o) arc asymptotic to negatively curved FLRW dust cosmology, as the latter 
follows from these asymptotic convergence forms. 

k q i — » and m q i ~ too. For models with a simultaneous big bang having these 
convergence forms we need to find the limit of Fh and F e as k q i — > with m q i = ttiq 
fixed: 

2 

lim Fh(\k q i\,m ) = lim F e (k qi ,m ) = — -. (118) 
|fc„i|->o k qi ^o 3y2TOo 

which fixes the constant bang time in terms of the level curve c(ij — t(a)) = 
2/ (3\/2mo). As a consequence, hyperbolic and elliptic models with t bb = 
and these initial value functions are asymptotic to a spatially flat FLRW dust 
cosmology, as the latter corresponds to these asymptotic convergence forms of 
m q i and k q i. 

mqi ~ > and k q i ~ fco- This asymptotic convergence only occurs in hyperbolic 
models (k q i < 0), since in elliptic models < k q i < 2m q i must hold. To examine 
this case in models with a simultaneous big bang, we evaluate the limit of Fh as 
iriqi with k q i — fco fixed: 



lim F e (k ,m ql ) 



(119) 



which fixes the constant bang time in terms of the level curve c(ij — t( )) = 
l/y/ko. Thus, hyperbolic models with t bb = t/ ) an d these initial value functions 
are asymptotic to the Milne Universe, which corresponds to these asymptotic 
convergence forms of m q i and k q i. 

k q i — > and m q i — > 0. In order to examine this asymptotic convergence for 
models with a simultaneous big bang, we explore the behavior of Fh and F e in 
the direction in the (m q i, k q i) plane given by the ray k q i — aQm q i, where a$ > 
is a constant («o < 2 for elliptic models). Hence, Fh and F e along this direction 
are given by 



F h (\k ql \) = 



ap Z h (l/a ) 



F e (k qi ) = 



«o Z e (l/a ) 

1 1/2 : 



(120) 



Thus, in the limit as m q i — > and k q i - 
limit as k q i tends to zero) we have Fh 



(which reduces in this direction to the 
■> oo and F e — > oo. In fact, both Fh and 
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F e only diverge as (rn q i,k q i) — > (0,0). All this means that c(U — t(o)) 00 
for any choice of m q i , k q i in which both tend to zero asymptotically. As a 
consequence, m q i and k q i cannot be fit to any level curve of these functions, and 
so a simultaneous big bang is incompatible with both k q i — > and m q i — > as 
r —> 00. Since this radial asymptotic behavior of initial value functions comprises 
models asymptotic to Minkowski in generalized Milne coordinates, self-similar 
dust solutions and Schwarzschild-Kruskal, models with a simultaneous big bang 
cannot be asymptotic to any of these spacetimes. This is not surprising since 
ct bb — > —00 holds in all these cases. 




Figure 1. Regularity domain of parabolic and hyperbolic models. The 

figure displays the (ct, r) plane with the locus of the initial curvature singularity 
(big bang) L = 0, marked by the thick curve [ct bb (r) , r] , together with the 
initial slice 3 T[<i] (dotted line). For all models asymptotic to a FLRW or Milne 
spacetimes (panel (a)), ci b b tends asymptotically to a constant, while models 
asymptotic to Minkowski in generalized Milne coordinates, self-similar solutions 
or Schwarzschild-Kruskal, we have ct bb — > —00 in this limit. Notice that space 
slices s T[t] are fully regular only for ct > B. 



12.2. Simultaneous maximal expansion. 



Let t m „ = t 



(m) 



denote the constant time value associated with "H„ = 0, then (91) 



implies that m q i and k q i are necessarily linked by the constraint 



"(t(m) ~ U) = F(k qi ,m qi ) 



_ rn q i [tt - Z e (k qi /m q i)} 

~ i, 3 / 2 : 

K qi 



(121) 



where Z e follows from ( 89 ) and we have expressed xi , yi in terms of m q i and k q i by 



means of (90a). This constraint is very similar to (1176), hence the compatibility 



with assumptions on radial convergence of m q j an d k qi follows from a qualitative 
study of the level curves of the function F in (121). As shown in figure 3, the level 
curves of this function do not allow for k q i — > to occur with m q i — > tuq > 0, 
since k q i — ¥ necessarily implies m q i —¥ 0. Therefore, models with in which ci max is 
constant must exhibit the same asymptotic states as the models examined in section 
11.5 (xj = k q i/m q i —> 2), in which ci max tends to a constant as r — > 00. In fact, the 
models of section 11.5 asymptotically converge to models with a simultaneous t mMt . 
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Figure 2. Regularity domain of open elliptic models. The figure displays 
the (ct, r) plane with the locus of the initial and collapse curvature singularities 
corresponding to L = 0, marked by the thick curves [ctbb(0, r ) an d [ct co ii{r), r], 
together with the initial slice 3 T[ti] (dotted line) and the maximal expansion time 
(H q = 0) marked by the thick dotted curve [rf max (r), r]. Panel (a) corresponds 
to models asymptotic to the spatially flat FLRW cosmology, so that ctbb tends 
asymptotically to a constant (while it tends to — oo in all other cases). Panel (c) 
corresponds to the case k q i/m q i — > 2 examined in section 11.5, for which ct max 
tends asymptotically to a constant. Panel (b) corresponds to all other cases. For 
the collapse time we have ct co n — > oo in all cases. Notice that in cases (a) and (b) 
the full asymptotic radial range is contained in the expanding phase (ct < ct max ). 
Also, the space slices 3 T[t] are fully regular only for B < ct < C. 




Figure 3. Level curves of F in p21| >. The fi gure displays the level curves in 
the (k q i,m q i) plane. We have F = in the line k q i = 2m gl , while F — > oo 
as k q i — > with m q i > 0, while all curves converge at (k q i,m q i) — > (0,0). 
As a consequence, only the asymptotic convergence k q i —> and m q i — > is 
possible, and so elliptic models with a simultaneous maximal expansion cannot 
be asymptotic to a FLRW cosmology. 
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Depending on the power law decay of m q i and k q i, these models are all asymptotic 
to Minkowski in generalized Milne coordinates, or to self-similar or Schwarzschild- 
Kruskal spacetimes. 



13. Summary, conclusion and final discussion. 

We have examined the asymptotic regime along the radial direction for regular 
parabolic, hyperbolic and open LTB models admitting a symmetry center and 
a covariant time slicing that defines a one-parameter family of 3-dimensional 
hypersurfaces, 3 T[t], orthogonal to the 4-velocity and marked by constant t. With 
the help of suitably defined covariant quasi-local scalars, m q , fc g , 1-L q and their relative 
fluctuations, this time slicing allows us to study LTB models under a well posed initial 
value framework in terms of a fiducial initial hypersurface 3 7~[ii] marked by t = tj. 
Since radial rays in the 3 T[i] are spacelike geodesies of the LTB metric whose affinc 
parameter is proper length, £, asymptotic conditions naturally correspond to the limit 
£ — > oo, which is related to the limit R — > oo along the radial rays (see section 5). 
However, £ = £{t, r) and R = R(t, r), and thus it is very difficult to use these invariant 
quantities to probe the behavior of covariant scalars in this asymptotic limit. Hence, 
we have provided in section 4 the appropriate conditions in which £ — ¥ oo corresponds 
to the limit r — > oo, so that a regular asymptotic regime follows by demanding that 
covariant scalars (local and quasi-local) remain bounded in this limit. The asymptotic 
limit along radial rays then follows from looking at the uniform convergence of initial 
value functions m q i{r), k qi (r) to analytic trial functions rh qi (r), k qi (r) as r — > oo, 



under the assumption of the radial coordinate gauge ( 36 ) , which fixes the radial 



coordinate as proportional to Ri = R(ti,r) (sections 6 and 7). 



13.1. Classification in terms of their asymptotic limits and states. 

By assuming asymptotic convergence forms rh q i, k ql for m q i, k q i, we examined sep- 
arately in sections 9-12 the radial asymptotic behavior of the metric functions and 
the covariant scalars, for parabolic, hyperbolic, open elliptic models and models with 
a simultaneous big bang or maximal expansion, making the distinction between the 
"asymptotic limit" as the strict limit as r — > oo, and the "asymptotic state" as the 
set of asymptotic expansions of all scalars around the asymptotic limit up to leading 
terms on rh q i(r), k q i(r) (see section 8). The summary of the results is as follows: 



Asymptotic Limit. All open LTB models can be classified in terms of two broad 
classes: asymptotic limit to FLRW or Minkowski. Considering A = to, H, k and 
q , 1-L qi k q , the "conventional" variables M, E and ct hh (bang time), as well as 



in 



ct,. aU , ci max (collapse and maximal expansion times, see (91 )), together with the metric 



functions L and F in (22), and using the radial coordinate gauge (36), we have: 



• Models whose asymptotic limit is a FLRW cosmology, including Milne Universe 
as a particular vacuum case (see figures la and 2a): 



L -> L(t), 
A q -> A q (t), 
M - too-RoT- 3 
ciooii — > oo, 



R -> Ro r L(t) -> oo, r -> 1, 
A^A q (t), S iA) ^0, 
E ~ ±fc -Ro r2 ' c *bb -» const., 
ci max — > oo (only elliptic models). 



(122a) 
(1226) 
(122c) 
(122 d) 
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• Models whose asymptotic limit is Minkowski, containing self-similar, 
Schwarzschild-Kruskal and generalized Milne asymptotic states (see figures lb 
and 2b-2c): 

L->1, R^R r ->oo, T -> 1, (123a) 
A,->0, A ~ A q (l + 5^) ->• 0, 5 {A) -> ^ (in general), (1236) 
M~m gi r 3 i^, E~k ql r 2 R 2 Q , ct bb -> -oo, (123c) 
c£ C oii — >• oo, ct maK — > oo (or — > const) (only elliptic models). (123d) 



where the special case ct max — > const, in ( 123a| has been discussed in sections 11.5 
and 12.2. 

Asymptotic States. It is useful to list these states for each kinematic class: 

Parabolic models (section 9) are radially asymptotic to 

• FLRW dust with zero spatial curvature (section 9.1) 

• Minkowski in generalized Milne coordinates (section 9.2) 

• Self-similar dust solution with zero spatial curvature (section 9.2) 

• Schwarzschild-Kuskal in Lemaitre coordinates, built by radial geodesies with zero 
binding energy (section 9.2) 

Hyperbolic models (section 10) are radially asymptotic to 

• A subclass of parabolic models (section 10.3), which are radially asymptotic to 

— FLRW dust with zero spatial curvature 

— Minkowski in generalized Milne coordinates 

• Milne Universe (section 10.4) 

• FLRW dust with negative spatial curvature (section 10.5) 

• Minkowski in generalized Milne coordinates (sections 10.4 and 10.5) 

• Self-similar dust solution with negative spatial curvature (section 10.5) 

• Schwarzschild-Kuskal in coordinates given by radial geodesies with positive 
binding energy (section 10.4) 

Open elliptic models (section 11) are radially asymptotic to 

• A subclass of parabolic models (section 11.3), which are radially asymptotic to 

— FLRW dust with zero spatial curvature 

— Minkowski in generalized Milne coordinates 

• Minkowski in generalized Milne coordinates (sections 11.4 and 11.5) 

• Self-similar dust solution with positive spatial curvature (sections 11.4 and 11.5) 

• Schwarzschild-Kuskal in Novikov coordinates, built by radial geodesies with 
negative binding energy (sections 11.4 and 11.5) 

Models with a simultaneous big bang (hyperbolic and open elliptic) (section 
12.1) are radially asymptotic to 

• FLRW dust with zero spatial curvature (as in sections 10.3, 10.5 and 11.3) 
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• FLRW dust with negative spatial curvature (as in section 10.5) 

• Milne Universe (as in section 10.4) 

Models with a simultaneous maximal expansion (elliptic) (section 12.2) are 
radially asymptotic to the same spacetimes as those of section 11.5. 

13.2. Cosmological considerations vs radial asymptotics. 

It is usual to consider any given LTB configuration as a model of a spherical 
inhomogeneity somehow immersed in a cosmic "background" . In order to discuss 
this issue in terms of the results obtained in this article, it is illustrative to relate the 
asymptotic forms of LTB scalars with observational parameters of a FLRW cosmology. 

If we denote by fi(f) the FLRW Omega parameter, a possible generalization for 
LTB models can be given in terms of quasi-local scalars: 

- m\ K 2 q 2m q - k q 2m ql - Lk ql ' 1 ' 

so that each kinematic class (parabolic, hyperbolic or elliptic) follows from the sign of 
fi - 1 

n ~ l = W q = 2m q k "~Lk qi - (125) 

While Cl is a covariant quantity (because m q and T-L q are covariant), it is not a quasi- 
local scalar: notice that we do not use the symbol VL q because £1 cannot be obtained 
from applying Q to the ratio 2m/H 2 . In fact, other expressions for the Omega and 
Hubble factors for LTB models have been suggested in the literature [39 , 40J (see [3U] 
for further discussion on this issue). However, it is evident that for all LTB models 
whose radial asymptotic limits and states are FLRW cosmologies we have as r — > oc 

n - n(t), (126) 

which clearly suggests that LTB models with this radial asymptotic behavior are 
suitable to model cosmological inhomogeneities asymptotically converging to a FLRW 
"background". This means that, irrespective of how we choose to generalize the 
observational parameters, they will tend asymptotically in the radial direction to 
the FLRW parameters. Notice that we have reached this conclusion by looking 
at the asymptotic states of LTB models defined all along the radial rays, without 
considering the rather artificial situation in which a FLRW radial asymptotic state 
of an LTB model is forced by matching to it to a FLRW cosmology at a finite fixed 
comoving boundary r = r^. Though such configurations can also be constructed 
with LTB models, leading to smooth and fully relativistic generalizations of the 
Newtonian "spherical collapse" or "top hat" models used in qualitative studies of 
structure formation (see [29] for discussion and references on this type of models and 
their relativistic generalizations). 

The cosmological interpretation of the radial asymptotics of LTB models with 
asymptotic limit to Minkowski is more complicated. To discuss these cases it is 
convenient to introduce a "local' equivalent of (l defined with the local scalars 
771 = np/3 and H — 0/3 as 

- _ K P _ 2777 (n-l)(l + 3(*))-(tfW) 2 
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where we eliminated % and k from (18 a I and used the constraint 

U 2 = 2m - k + S 2 = 2m - k + {U q 5 [H) ) 2 , (128) 

which follows from the definitions of G and X in ^ and Q. For models with an 
asymptotic FLRW limit, -> and <jW -> hold, so ^i oc ->■ 0, and the FLRW 
fl(t) results. However, in models with a Minkowski limit do not tend to 

zero, resulting in a different form for Cli oc . 

In parabolic models with asymptotic li mit t o Minkwoski we have fi = 1 
everywhere, while 8^ — > 6+ = 5^/2 (from (506)), hence 

n loc ^i-i--^ m \ <i, (129) 

Since we have ft — >• 1 for hyperbolic with a matter dominated Minkowski limit (section 
10.4) and for all open elliptic models with a Minkowski limit (which is also matter 
dominated), we obtain for these cases the same limit as in ( |129[ ). For all hyperbolic 
models with a vacuum dominated Minkowski limit we have £~2 — >■ 0, hence 

2<$(?i) - J( fc ) 

(1+J(«))2 ■ (130) 

Therefore, considering ((75rf|-((75e| , (|78c|)-(|78^) and (l87c|-((87d|), we have 5<*0 -> 



25^\ and so (l\ oc — > holds for models asymptotic to Milne, generalized Milne and 
self-similar solut ions, b ut in models asymptotic to Schwarzschild-Kruskal we have 



(129) 



<j( fe ) _^ ( see (8Xd)-(81e|), so for these models Cl\ oc has an asymptotic form as in 



Since the asymptotic forms for either Cl or Cl\ oc are independent of t (being either 
zero or a constant < 1), a cosmological interpretation for the radial asymptotics of 
models with a Minkowski limit cannot be given in terms of a smooth convergence 
or transition to a FLRW model in the radial direction. Instead of a cosmic FLRW 
background, the external realm of these LTB models could be a large void region. 
A similar and appealing interpretation of radial asymptotics to Minkowski follows 
naturally from the notion of "finite infinite" (fi), introduced by Ellis [37 (see also 
comments in [38 ) to represent a near Minkowskian timelike envelope around bound 
structures, which defines a scale in which these structures can be studied as almost 
asymptotically flat systems without considering the effect of cosmic expansion. Such 
a scale would be intermediate between characteristic lengths of bound structures and 
a cosmic scale where expansion cannot be ignored. Evidently, LTB models with 
asymptotic limit to Minkowski approximate spherically symmetric realizations of this 
cosmic structure in such a near Minkowski envelope. 

The results of this article are useful for probing supernovae and CMB observations 
by means of LTB models. They are essential for the study of radial profiles of covariant 
scalars, namely: the "clump" or "void" profiles and the possibility of profile inversions 
as the models evolve in time. This study has been undertaken in a separate article 
[41] . The radial asymptotic properties of the models have also theoretical and practical 
consequences in the application of Buchert's scalar averaging formalism |23] to LTB 
models, as vacuum dominated LTB models are the LTB configurations most likely to 
yield an "effective" acceleration that mimics dark energy emerging from back-reaction 
terms in the context this formalism. These connections have been already remarked 
[211 122j and are currently under further investigation. 
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Appendix A. Analytic solutions in the conventional variables. 

The solutions of the Friedman-like equation ^ in terms of the conventinal free 
functions M, E, ct bb are: 

• Parabolic models or regions: E = 0. 

c(t-t bb ) = ^r, 3 , R=(2M)^ V 2 , (A.l) 

• Hyperbolic models or regions: E > 0. 

M M 
R = — (cosh n - 1) , c(t - t bb ) = (sinh r, - rj) (A.2) 

Z?3/2 Ej R 

— c(U-t hb )=Z h (U), U=—, (A.3) 

• Elliptic models or regions: E < 0. 

i? = (1 - cos 77) , c(i - t bb ) = (»7 - sin rj) , (A.4) 

|£| 3/2 , Z e {U) expanding phase R > _ |^|-R 

M C[ t hb) ~ { 2TT-Z e (U) collapsing phase R < ' ~ M ' 

(A.5) 



where Z/i and Z e are respectively given by ( 52 ) and ( 89 1 . The solutions presented 



and (88). 



above are given in terms of the new variables by (41 1, (51 ) 

Another particular solution of ([2| follows by assuming M = and E > (but 
otherwise arbitrary), leading to the case "[s2]" in [S] 

R=Ec(t-t bb ), (A.6) 

which are locally Minkowski solutions (sections of Minkowski spacetime in coordinates 
that generalize Milne's Universe). 



Appendix B. Prescribing local initial value functions. 



In section 7 we assumed an asymptotic convergence for for A qi and obtained the forms 

for Ai and 5^ . We examine here the opposite situation. If we assume that A4 ~ A4 , 
then for r > y and with the help from (|9| and ( 36 ) we have 



A — 



" ry 




/ A; t x 2 Ax + 


/ ylia; 2 da: 


7o 





(B.l) 



A^a: da: ~ 

while 5^ follows directly from |llj: 

{A) Ai (r)r 3 - 3 A z x 2 dx - 3 f y A lX 2 dx 

1 3 J y A x x 2 dx + 3 f y A t x 2 dx ' 

A comparison of (38a)- ([386| and ( |B.1| ( B.2[ ) shows that it is much easier to obtain the 
asymptotic limit of J- from A q i than from A4. Prescribing A q i is also more useful 
in practice, since the initial value functions in the analytic solutions given in sections 
9-11 are m q i, k q i. 



(B.2) 
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We now examine the asymptotic forms for m qi and d^ 71 " 1 that follow given a 



prescribed asymptotic form for mj. Consider a power law decay m, 



m r 



(B.ll and (B.2) yield 



' (jl 



m 



Am) 



(3- a)r c 
where: I(y) 



3 

m 



-ar~ a - (3- a)I{y)r~ 3 
{3 - a)I{y)r- 3 +3r- a 

3 1 — 



then 



(B.3) 



rriiX 2 dx - 



-y 



iq o — a 

Notice that if a < 3 then r~ a is the dominant term, whereas if a > 3 then r~ 3 is 
dominant. Hence, we obtain m qi oc r~ a for a < 3 and for a logarithmic decay, while 
m q i oc r -3 results if we assume that a > 3 or an exponential decay. However, from 
Lemma 3, both rrii and m q i tend to zero. By assuming r>i;we obtain the following 

asymptotic limit for 5^ 



r(m) 



-a/3, 
-1 



if 
if 



a < 3 
a > 3 



(B.4) 



(m) 



-1, 



It is straightforward to show that if decays to zero exponentially, then 5, 
while a slow logarithmic decay yields 5^ ' — > 0. As opposed to prescribing an 
asymptotic decay for m q i — > as in (39a)-(39&l, all assumptions on the decay of 
necessarily yield m,; and m q i positive in the full range r > y. This is so 



m 



because m q i > > holds in this range (as m! i and m' qi are negative, see (10a)) 



Appendix C. Particular cases of open LTB models. 



Given the radial coordinate gauge (36), all open LTB models can be univocally 
characterized by the initial value functions m q i, k qi . Hence, spacetimes that emerge 
as particular cases of these models emerge by specific specialization of these functions. 



The line elements of these spacetimes follow directly from ( 22 ) by recalling the 



particular cases of these functions accordingly, and define particular case spacetimes 
which have served as references to the possible asymptotic states of all open LTB 
models in the limit r — >• oo. 



Appendix C.l. FLRW dust spacetimes 

Homogeneous and isotropic dust FLRW cosmologies follow as particular cases of LTB 
models by setting 



m > and fc„ 



where m , k are constants. From (21), (42), (54) and (91), we have 

„3 



M 



m r 



E 



-fc r 2 , 



constants. 



(C.l) 



(C.2) 



Evidently, the spatially flat case corresponds to fco = 0> while positive/negative kg are 
FLRW models with positive/negative spatial curvature. It is evident from (41) that 
fl(Trj ) yields for k = 

2/3 



L(t) 



1 



V2m c(t - U) 



fco = 0, 



(C.3) 
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while for fco ^ the lunction L — L(t) follows from the implicit solutions (511 and 
under the specialization (C.l): 

|*o | |fco| 3 / 2 

> 2/o = , 

mo m 

with Z 



<P(t) = 2/o c(t - U) + Z(x ) = Z(x L), 



x 



(C.4) 



An) 



Zh or Z = Z e given by ( |52| and ( 89 ) for hyperbolic (negative curvature) and 
elliptic (p ositiv e curvature) FLRW models. For all particular cases (C.l), equations 
fllf and |l8a| imply 6$ 
Hi - 

S (m) = §(k) = 0; and SOj g(ftr= o, m = 



0, so that rrii 



k„i and 



U ql hold. Also, (15), (43) and (55) imply T = 1, and thus (19ol and (196) yield 



i q , k = k q and H = T-L q follow from (17a 



(176) and (20). As a consequence, (22) becomes a FLRW line element and all scalars 



are either constants or functions of t only. 



Appendix C.2. Minkowskian particular cases. 

Considering the Riemann tensor (see Appendix A2 of [30 ) , any particular LTB model 
in which m q — m = is locally equivalent to Minkowski spacetime given in curvilinear 
coordinates that generalize the well known Milne Universe. These are the locally 
Minkowskian solutions (A. 6), corresponding to the case "[s2]" in the classification of 



[8], and are given by the specialization: 



and 



k qi < 



but otherwise arbitrary. 



(C.5) 



so that m q = m = = 0, but 5^ ^ 0, 5^ ^ 0. The metric functions L and T 



definition of T in (12) 



follow from substituting (C.5) into the equation for L in (15) and substituting in the 



l + \k qi \^ 2 c (t-U) 7 



+ M 1/2 (l + l^) c(t-t,) 



1 



(C.6) 



where we used (11 1 and the coordinate gauge (36). Evidently, wc have M = 0, while E 



is arbitrary as in (A.6). The scalars associated with spatial curvature and expansion 
correspond to hypersurfaces 3 T and a 4-velocity associated with test observers. In 
general, these coordinates do not cover the full Minkowski spacetime, and thus L = 0, 
marked by ct — ct hb (r) — cti — l/l/e^l 1 / 2 , corresponds to a caustic surface for this 
congruence, not to a curvature singularity. Milne's Universe, which can also be 



considered as a vacuum FLRW cosmology, is the particular case of (C.5) and (C.6) 
given by 

m q i = and \k q i\ — ko = constant, (C.7) 

L = L(t) = l + y/ktoc(t-ti), (C.8) 

M = 0, E = T 2 - 1 = |fc |r 2 , t bb = constant. (C.9) 

in this case vanish, then T — 1 and all the 5^ also vanish. Also, we 
n = for all t, though k q = k < and % q = % are functions of time. In 



Since 5^ 
have m q = 

the general case (|C.5| k q = k < and H q = H are also functions of r. 



Appendix C.3. Schwarzschild-Kruskal spacetime. 

The Schwarzschild-Kruskal spacetime in comoving coordinates constructed by radial 
geodesies (Lemaitre and Novikov coordinates, see page 332 of [3]) follows by setting 

m q i — mo r~ 3 and k q i arbitrary, (C.10) 
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Hence, (21) and (361 imply 

M = M = too #o. 
while {!]} an d |l9a| yield = 6^ 



E = F 2 - I arbitrary, (C.ll) 

— 1, so that 2m g = 2Mq/R 3 and to,; = m = 
follows from (|18a| and (|17a[). It is straightforward to prove that the Riemann tensor 



(see Appendix A2 of |30]) with to = and m q ^ leads to a vanishing Ricci tensor, 
and thus we can identify this case as a vacuum Schwarzschild-Kruskal solution where 
2Mo is the Schwarzschild radius. The variables associated with spatial curvature 
(k, k q and 5^) and the expansion (U, n q and 5^) remain arbitrary because a 
spatial curvature and an expansion can always be associated to a congruence of radial 
geodesic test observers with arbitrary binding energy. In fact, by setting M = M in 
([2]) we obtain the equation for radial timelike geodesies in a Schwarzschild spacetime 
if we identify J- 2 — 1 with the binding energy of the comoving test observers. 

The form of L depends on the choice of k q i. If we choose geodesic observers with 
zero binding energy k q i = 0, then L and t bb follow from (41) and (42) with m q i given 



by (C.10I 



L(t,r) 



1 + -V2m r- 3/2 c(t-U) 



2/3 



Cti 



2r 3/2 

3-^/27710 



(C.12) 



If we choose k q i ^ 0, then L follows from the implicit solutions (51) and (88) under 



the specialization (C.10) 



4>(t, r) = yi c(t - U) + Z(xi) = Z(xi L), 



k q i T 3 

to 



7,3/2 3 

K qi r 



to 



(C.13) 



where Z is either one of the functions ( 52 ) or ( 89 ). Notice that t 

(and t, 



—00 as r 



(C.12), whereas the form of t b 



As shown in sections 10 and 11, if k oi r 



qt 



t coll if k qi > 0) depend on the choice of k qi 
■> 0, then L tends asymptotically to (C.12) 



We also remark that the Schwarzschild-Kruskal spacetimer lacks a regular symmetry 
center, hence r = is not a special worldline. Also, the comoving coordinates in the 
cases k q i < do not cover the full Schwarzschild-Kruskal manifold. 

Appendix C.4- Self-similar solutions. 

The self-similar LTB solutions have been popular models to study gravitational 
collapse [TU]. They have been classified in |36j . In terms of the parameters used 
in this article these solutions follow as the particular case 

m q i = mor~ 2 and k q i = k^r~ 2 1 (C.14) 



where toq > and ko are constants. The form of L follows from (41) and from the 



implicit solutions (|51|) and 



1 2/3 



under the specialization (C.14) 
o=0 



$(t, r ) = Vq C + Z(xq) = Z(x L), 



x 



k 
m 



7-3/2 
^0 

toq : 



(C.15) 



(C.16) 



where Z = Zh or Z = Z e given by (|52|) and ((89 1 for hyperbolic and elliptic models, 
C='^- — ^. (C.17) 



and we can identify the self similar variable 
c(t - U) 
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The functions M, E and t b 
M = m Rl r, E = T 2 

k = 0, 



follow from applying (|C.14|) to (|2l), (|42]) and (|54| 

- 1 



(C.18) 
(C.19) 



-oo as r 



oo (see section 13). 
-2/3. The forms for 



2r 

3 V2m ' 2/o 
As in the Schwarzsch ild-Kr uskal case, we have t bb - 
As a consequence of (C.14) we have (in general) 8[ 

m q , 1-L q , k q , r, 5( m \ 6^ k> , and the remaining scalars follow by applying (C.14) to 
the appropriate equations. It is straightforward to verify that the metric (22 1 takes the 
appropriate self-similar form, admitting the nomothetic vector z a = c(t — ti)SQ + n5". 
Also, all dimensionless variables (i.e. of the form A/Ai, as well as L, V and W q c(i— 
depend only on the self-similar variable. 



References 



[1] Lemaitre G 1933 Ann. Soc. Sci. Brux. A 53 51. See reprint in Lemaitre G 1997 Gen. Rel. Grav. 
29 5; Tolman R C 1934 Proc. Natl Acad. Sci. 20 169; Bondi H 1947 Mon. Not. R. Astron. 
Soc. 107 410. 

[2] Krasinski A, Inhomogeneous Cosmological Models, Cambridge University Press, 1998. 

[3] Plebanski J and Krasinski A, An Introduction to General Relativity and Cosmology , Cambridge 

University Press, 2006. 
[4] Krasinski A and Hellaby C 2002 Phys Rev D 65 023501 
[5] Krasinski A and Hellaby C 2004 Phys Rev D 69 023502 
[6] Krasinski A and Hellaby C 2004 Phys Rev D 69 043502 
[7] Hellaby C and Krasinski A 2006 Phys Rev D 73 023518 

[8] Matravers D R and Humphreys N P 2001 Gen. Rel. Grav. 33 53152; Humphreys N P, Maartens 

R and Matravers D R 1998 Regular spherical dust spacetimes Preprint gr-qc/9804023vl. 
[9] Bolejko K, Krasinski A and Hellaby C 2005 MNRAS 362 213-228 
[10] Eardley D M 1974 Commun Math Phys 37 287; Eardley D M and Smarr L 1979 Phys Rev D 

19 2239; Dyer C C 1979 MNRAS 189 189; Waugh B and Lake K 1988 Phys Rev D 38 1315; 

Waugh B and Lake K 1989 Phys Rev D 40 2137; Lemos J P S 1991 Phys Lett A 158 279 
[11] Joshi P S and Dwivedi I H 1993 Phys Rev D 47 5357; Joshi P S and Singh T P 1995 Phys Rev 

D 51 6778; Dwivedi I H and Joshi P S 1997 Class. Quant. Grav. 47 5357 
[12] Vaz C, Witten L and Singh T P 2001 Phys Rev D 63 104020; Kiefer C, Mueller-Hill, Vaz C 2006 

Phys Rev D 73 044025; Bojowald M, Harada T and Tibrewala R 2008 Phys Rev D 78 064057 
[13] Pascual-Sanchez J F 1999 Mod. Phys. Lett. A 14 1539; Sugiura N K and Harada T 1999 Phys 

Rev D 60 103508; Celerier M N 2000 Astron. Astrophys. 353 63; Tomita K 2001 MNRAS 

326 287; Iguchi H, Nakamura T and Nakao K 2002 Prog. Theor. Phys. 108 809; Schw arz D 

J 2002 Accelerated expansion without dark energy Preprint arXiv : astro-ph/0209584v2; 
[14] Apostolopoulos P et al 2006 JCAP P06 009; Kai T, Kozaki H, Nakao K, Na mbu Y and Yoo C 

M 2007 Prog. Theor. Phys. 117 229-240 {Preprint |arXiv:gr-qc/0605120| l; Mattsson T and 

Ronkainen M 2008 JCAP 0802 004 {Prep rint arXiv : astro- ph/0708 . 3673v2) ; Bolejko K and 

Andcrsson L 2008 JCAP 10 003 ( Preprint [arXi v : 0807 . 3577) 
[15] K olb E W, Matarrese S, Notari A and Riotto A 2005 Phys Rev D 71 023524 {Preprint 

|arXlv : hep = ph /0409038v2) ; Marra V, Kolb E W and Matarrese S 2008 Phys Rev D 77 023003; 

Marra V, Kolb E W, Matarrese S and Riotto A 2007 Phys Rev D 76 123004. 
[16] Garcfa-Bellido J and Troels H 2008 JCAP 0804:003 Preprint gr-qc/0802 . 1523v3 [astro-ph] 
[17] Moffat J W 2006 J. Cosmol. Astropart. Phys. JCAP(2006)001; Alnes H, Amazguioui M and 

Gron O 2006 Phys Rev D 73 083519; Alnes H and Amazguioui M 2006 Phys Rev D 74 

103520; Alnes H and Amazguioui M 2006 Phys Rev D 75 023506 
[18] Rasanen S 2006 Class . Quant. Grav. 23 1823-1 835; Enqvist K and Mattsson T 2007 JCAP 

0702 019 (Prepr mi|arXiv:astro I ph/0609120/4); Enqvist K 2008 Gen. Rel. Grav. 40 451- 

466 (Preprint arXiv : 0709 . 2044 ) 
[19] C huang C H, Gu J A and Hwang W Y P 2005 Class. Quant.Grav.,25, 175001 Preprint 

astro-ph/0512651 
[20] Paranjape A and Singh T P 2006 Class. Quant.Grav., 23, 69556969 

[21] Sussman R A 2008 On spatial volume averaging in Lematre-Tolman-Bondi dust models. Part 
I: back reaction, spacial curvature and binding energy Preprint arXiv: 0807 . 1145 



Radial asymptotics of LTB dust models. 



44 



[22] Sussman R A 2009 Quasi-local variables and scalar averaging in LTB dust models Preprint 

I arXiv: 0912. 40741 

[23] Buchert T 2000 Gen. Rel. Grav 9 306-321 (Preprro* |arXiv:gr- qc/0001056]rl); Buchert T 2008 

Gen. Rel. Grav. 40, 467 
[24] Celerier M N 2007 New Advances in Physics 1 29 (Preprmt |arXlv : ast"ro 1 ph/0702416 ) 
[25] Wainwright J and Andrews S 2009 Class. Quant.Grav., 26, 085017 
[26] Sussman R A and Garcfa-Trujillo L 2002 Class. Quant.Grav. 19 2897-2925. 

[27] Sussman R A Quasi-local variables and inhomogeneous cosmological sources with spherical 

symmetry 2008 AIP Conf.Proc. 1083 228-235 Preprint |arXiv:0810. 1120| 
[28] Sussman R A 2008 Class Quantum Grav. 25 015012 Preprint arXiv : grqc/0709 . 1005 
[29] Sussman R A 2009 Phys Rev D 79 025009. 

[30] Sussman R A 2010 A new approach for doing theoretical and numeric work with Lemaitre- 

Tolman— Bondi dust models Preprint arXiv: 1001 . 0904vl 
[31] Ellis G F R and Bruni M 1989 Phys Rev D 40 1804 

[32] El lis G F R and v an Elst H 1998 Cosmological Models (Cargese Lectures 1998) Preprint arXiv 
gr-qc/9812046 v4 

[33] va n Elst H and Ellis G F R 1996 Class Quantum Grav 13 1099-1128 (Preprint 
arXiv:gr-qc/9510044) 

[34] Hayward S A 1996 Phys Rev D 53 1938 (Preprint Ar Xiv |gr-q c/9408002); Hayward S A 1998 
Class Quantum Grav 15 31473162 (Preprint ArXiv gr-qc/9710089v2) 

[35] Hellaby C and Lake K 1985 Astrophys J. 290 381 

[36] Carr B J 2000 Phys Rev D 62 044062 

[37] Ellis G F R, in Bertotti B de Felice F and Pacolini A (eds), General Relativity and Gravitation 

(Reidel, Dordrecht, 1984) pp 215-288; 
[38] Wiltshire D 2007 New J. Physics 9 377 (Preprint ArXiv |gr-qc/ 0702082v4) 
[39] Moffat J W and Tatarski D C 1995, Astroph J. 453 17. 

[40] Humphreys N, Maartens R and Matravers D 1997 Astroph J. 477 47. See also Ehlers J 1993, 

Gen. Rel. Grav. 25 1225 (translation of original 1961 article). 
[41] Sussman R A 2010 Radial profiles of covariant scalars in Lemaitre-Tolman— Bondi dust models 

(in preparation). 



